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Abstract 

We present a novel analysis of the irN scattering amplitude in covariant baryon chiral per- 
turbation theory up to 0(p 3 ) within the extended-on-mass-shell renormalization scheme 
and including the A(1232) explicitly in the 5-counting. We take the hadronic phase shifts 
provided by partial wave analyses as basic experimental information to fix the low-energy 
constants. Subsequently, we study in detail the various observables and low-energy theo- 
rems related to the ttN scattering amplitude. In particular, we discuss the results and chiral 
expansion of the phase shifts, the threshold coefficients, the Goldberger-Treiman relation, 
the pion-nucleon sigma term and the extrapolation onto the subthreshold region. The chi- 
ral representation of the amplitude in the theory with the A presents a good convergence 
from very low energies in the subthreshold region up to energies well above threshold, lead- 
ing also to a phenomenological description perfectly consistent with the one reported by 
the respective partial wave analyses and independent determinations. We conclude that a 
model-independent and systematic framework to analyze 7r TV-scattering observables using 
directly experimental data shall be possible in covariant baryon chiral perturbation theory. 



1. Introduction 

Pion-nucleon (jrN) scattering is a hadronic reaction that gives access to some of the most 
prominent and fundamental questions related to the strong interactions [lj]. At low energies, 
it allows to test the dynamical constraints imposed by the chiral symmetry of QCD in one of 
the simplest processes involving a nucleon ^ . Also, understanding irN scattering is essential 
for a first-principle approach to the nuclear structure,given that the long-range part of the 
iViV interactions is given by the exchange of pions [3(. Experimental data on differential 
cross sections and polarization observables have been collected in the last 50 years, and more 
intensively in the last decade thanks to fully dedicated experiments run in meson factories. 
The usual way to organize the experimental information is by means of energy-dependent 
parameterizations of the scattering amplitude projected in partial waves, fitted to the data 
and supplemented with unitarity and analyticity constraints. These partial wave analyses 
(PWAs) provide an accurate representation of the data included in the parameterizations, 
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which can be used to extract values of scattering parameters and strong coupling constants 
or to identify the effect of resonances in the different isospin- angular momentum channels]]] 
In spite of the very long and sustained effort in studying the ttN scattering amplitude, 
there are fundamental questions concerning this process that have not been satisfactorily 
answered yet. For instance, it remains unclear what is the exact value of the pion-nucleon 
coupling constant g^N ° r of the pion-nucleon sigma term, a^N- Besides that, it is impor- 
tant to determine accurately the threshold parameters. The scattering lengths ought to be 
compared with those obtained from the analysis of the accurately measured IS* level shift in 
pionic hydrogen and deuterium jsj . An alternative source of phenomenological information 
on the scattering parameters is given by the analyses of the NN interaction 0|. On the 
other hand, it remains a challenge for the theory to understand hadronic processes directly 
from the parameters and dynamics of the underlying QCD. Important progress in the com- 
putation of hadron-hadron scattering lengths in unquenched LQCD has been reported jsj], 
although those of the nN system are still computationally prohibitive. Nevertheless, results 
were obtained in a pioneering quenched calculation ji| and, more recently, in an unquenched 



one 



10J. The LQCD simulations are often run with quark masses heavier than the physical 
ones and their results require a careful chiral extrapolation to the physical point. 

Chiral perturbation theory (ChPT), as the effective field theory of QCD at low ener- 



gies [111 |12L 1 131. 1 141] . is a suitable framework to build a model-independent representation of 



the nN scattering amplitude and to tackle some of these problems (see Ref. (15| for compre- 
hensive reviews). An interesting feature of ChPT is that, regardless of the specific values of 
its parameters or low-energy constants (LECs), it inherits the chiral Ward-Takashi identities 
of QCD among different Green tactions Q . Some of these identities at leading order in the 
chiral expansion were obtained using PCAC and Current- Algebra methods in the fifties and 
sixties, conforming what has been know as low-energy theorems since then [l6j]. Remarkable 
examples of these theorems in the ttN system are the Goldberger-Treiman (GT) relation 



between g nN and g& [17] and the one at the Cheng-Dashen (CD) point between a nN and 



the scattering amplitude [18|. Other interesting examples are the Weinberg predictions for 



the scattering lengths [19] or the Adler condition for the nN scattering amplitude 20 . 

ChPT allows one to investigate up to which extent the low-energy theorems apply since 
it provides a method to compute systematically the higher-orders in the chiral expansions of 
the correlation functions entering corresponding Ward-Takahashi identities. A first attempt 
to apply baryon ChPT (BChPT) to elastic irN scattering was undertaken by Gasser et al. 
in Ref. [ij], where an off-shell amplitude was obtained up to 0(p 3 ) in a manifestly Lorentz 



covariant formalism (for reviews on BChPT see Refs. [2l|, |22|). However, it was shown that 

the presence of the nucleon mass as a new large scale in the chiral limit invalidated naive 

power counting arguments in the baryon sector. This problem can be solved in the heavy- 

baryon formalism (HBChPT) jiij, in which one recovers a neat power-counting scheme at 

the cost of manifestly Lorentz covariance. Calculations of the irN scattering amplitude up 
,3"\ lo/ili OT1 J 



to 0(p ) [24] and 0(p ) [25] accuracies have been performed in the HB formalism by Fettes 



^ee the classical treatise of Hohler [l| for an exhaustive review on the nN scattering amplitude. For 
updated descriptions of the current experimental situation and data see Refs. 0, 01- 
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et al, showing a good description of the S- and P-wave phase shifts of different PWAs at 
low energies. 

It was latter shown that the HB approach is not well suited for studying some of the 



low-energy theorems involving the ttN scattering amplitude |26|. The problem is that the 



non-relativistic expansion implemented in HBChPT alters the analytical structure of the 
baryon propagator such that the chiral expansion of some Green functions does not converge 



in certain parts of the low-energy region [26|, |27J. This problem shows up in the analytic 



continuation of the nN scattering amplitude onto the subthreshold region [25| or in the 
behavior of the form factors close to the two-pion threshold, t = 21J. Besides that, it 

has been shown that the non-relativistic expansion may have a problematic convergence in 
some other cases 

Based on the ideas previously discussed in Ref. {Hj], Becher et al. proposed the in- 
frared scheme (IRChPT) [26| as a solution to the problems of the HB formalism. The IR 
scheme is a manifestly Lorentz covariant approach to BChPT that preserves the HB power 
counting at the same as it resumms the kinetic terms of the positive-energy part of the 
baryon propagators, curing the analyticity problems of the HB approach. The 7rN scatter- 



ing amplitude has been also calculated at 0(p ) [32|, [33J and 0(p ) [34J accuracies in the 
IR scheme. At 0(p 4 ) the amplitude rapidly converges in the proximity of the CD point 
so one can investigate meticulously the corresponding low-energy theorems. However, and 
similarly to HBChPT j25|, the convergence at 0(p A ) in IR is such that one fails to connect 
the subthreshold and threshold regions. 

A serious drawback of the IR method is that, in curing the problems of the HB expansion, 



it runs into its own ones with the analytic properties of the loop integrals [351 ] . This is related 



to the fact that the IR resummation of kinetic terms performed on the HB propagators 
completely omits the inclusion of negative-energy pole or anti-nucleon contribution, violating 
charge conjugation symmetry and, therefore, causality 36(. The most striking consequence 



of this is the appearance of unphysical cuts. Despite lying outside the range of applicability 
of ChPT, these cuts can have sizable contributions to the Green functions at low-energies, 
disrupting the convergence of the respective chiral expansions. This has been indeed shown 
in different applications of ChPT in the baryon sector such as the chiral extrapolation of the 
nucleon magnetic moment 28|, the SU(3)p breaking of the baryon magnetic moments 
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or the unitarized description of the ttN scattering amplitude [33 



A completely different difficulty in the baryon sector of ChPT is related to the A(1232)- 
resonance. Its contributions can be important at very low energies since this resonance is very 
close in mass to the nucleon. In the conventional chiral expansion, these effects are accounted 
for by the LECs but the radius of convergence, in this case, becomes drastically reduced. 
This problem is prominent in ttN scattering as the threshold for this process is at a center- 
of-mass (CM) energy 5W = M n away from the point around which the chiral expansion is 
performed. An improved convergence of the chiral series can be obtained including the A 



resonance as an explicit degree of freedom [37|, |38|, ]39j . In this case, one introduces a power 
counting for the new scale e = -mjy and computes the A contributions accordingly. 



The A corrections to the irN scattering amplitude have been calculated in the HB [40 
and IR 32] schemes up to 0(e 3 ) within the small-scale expansion (SSE), which considers 
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e ~ 0{p) [38j . In case of the HB calculation, the inclusion of the A increases the range of 
energies described as compared with the A-less case at 0(p 3 ), although the values of the 
LECs strongly depend on the fitted data, precluding a clear discussion on the extracted 
values of the observables related to nN scattering (e.g. a n pf). In the IR calculation, the A 



corrections worsen severely the description of the different PWA phase-shifts [32 



In this paper we present a novel chiral representation of the ttN scattering amplitude 
with two main differences as compared with previous work. In the first place, we use 
Lorentz covariant BChPT with a consistent power counting obtained via the extended-on- 



mass-shell (EOMS) renormalization scheme 4l|, |42|. The main advantage of EOMS over 



the HB and IR schemes is that it preserves the proper analytic structure of the Green 
functions j35|. Secondly, we explicitly include the A(1232) taking into account that, below 
the resonance region, the diagrams with the A are suppressed in comparison with those with 



the nucleon 43, 39 



The paper is organized as follows: We first present in Sec. |2]a detailed description of the 
formalism employed, which comprises the representations of the nN scattering amplitude 
and the ChPT formalism. Special emphasis is put on explaining the EOMS scheme (Sec l2.3p 
and the inclusion of A(1232) degrees of freedom (Sec. 12.41) . A brief discussion of isospin 
breaking effects is also presented here (Sec. I2.5p . In Sec. [3], we explain in some detail the 
calculation of the scattering amplitude and the application of the EOMS renormalization 
scheme. The reader is addressed to the Appendices for the complete results and some 
technical details. Sec. H] focuses on the determination of the LECs of BChPT without 
(^-ChPT) and with (A-ChPT) the contributions of the A (1232) as an explicit degree of 
freedom and using various phase-shifts sets provided by different PWAs. The convergence of 
the chiral expansion of the amplitude in either case is studied in Sees. 14.1.11 and 14.2.11 Once 
the LECs are determined, Sec.[5]is dedicated to discuss the results and the chiral structure of 
the different 7riV-scattering observables. The analysis of the threshold parameters, including 
a comparison with LQCD results, the GT relation and the pion-nucleon sigma term are then 
presented in Sees. 15. 1[ 15.21 and 15.31 respectively. Finally, we study the important issue of the 
connection between the threshold and subthreshold regions of the ttN scattering amplitude 
in Sec. |6j We close the paper with a summary of our findings and results together with a 
brief outlook. 



2. Formalism 

2.1. The scattering amplitude 

We consider the process n a (q)N(p, a; a) — > ir a ' (q')N(p' , a'; a'). Here a and a' denote the 
Cartesian coordinates in the isospin space of the initial and final pions with four-momentum 
q and q', respectively. Regarding the nucleons, cr(a') and a(a') correspond to the third- 
components of spin and isospin of the initial (final) states, in order. Since we calculate the 
on-shell amplitude, the usual Mandelstam variables s = (p + q) 2 = (p' + q') 2 , t = (q — q') 2 = 
(p — p') 2 and u = (p — q') 2 = (p' — q) 2 , fulfill s + t + u— 2M 2 + 2m 2 N , with m N and M n the 
nucleon and pion mass, respectively. Exact isospin symmetry is assumed in the following, so 
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Figure 1: Diagrammatic representation of the Lorentz- and isospin- invariant amplitudes. In this 
figure, s and t correspond to the Mandelstam variables, while P and P' correspond to the momen- 
tum of the incoming and outgoing nucleon, respectively. On the other hand, q (q') corresponds to 
the momentum of the incoming (outgoing) pion with Cartesian isospin index a (a 1 ). 



it is convenient to consider Lorentz- and isospin-invariant amplitudes. We then decompose 
the scattering amplitude as [l| 



T„n> — 5„' n T^ 



-[r a ,T a ,]T- 



T 



± 



u(p', a' 



u{p, a) 



where the Pauli matrices are indicated by r c and the definitions of the indices and momenta 
are shown in Fig. [TJ 

However, the decomposition (Q]) is not transparent from the power counting point of 
view, since the leading order contributions of and cancel. Therefore, it is convenient 
to write the scattering amplitude in terms of the energy variable v = J^-H that scales as 
v ~ 0(p), and the amplitudes D = A + vB and B. In this representation, the scattering 
amplitude takes the form 



u(p',a') 



D z 



1 



4m 



N 



u(p,a). 



(2) 



Furthermore, one can construct the representation in partial waves of the amplitude, which 
is obtained by projecting Eqs. ([!]) onto states with definite total angular momentum J and 



2 In fact, the chiral expansion can be understood as an expansion around the soft point, in which v 
0, t = and M n = 0. 
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orbital angular momentum £, 



TjAa! , a': a, a) = — = — 

y/4n(2e+l)(0<T<r\e±J) 

x^J dp' (7r(-p'; a')N(p', a'; a')\T\*(-p; a)N(p, a; a))(ma'a\ih)Yr(pT- 

m,cr' 

(3) 

This representation is often used in nN scattering phenomenology and, in fact, global 
analyses of the experimental database are often presented in terms of partial-wave phase 
shifts BHHS. In the |Appendix A] we specify the conventions used in this paper for the 



partial-wave representation of the scattering amplitude and the formulas used to extract the 
different phase shifts from the chiral representation. 

2.2. Chiral perturbation theory and chiral Lagrangians 

In ChPT one develops a power counting to express the Green functions of QCD as an 
expansion in the light quark mass and soft external momentum (generically denoted as 
p), which are small compared with the chiral symmetry breaking scale, A x ~ 4:irf n ~ 1 
GeV. These are called chiral expansions and they contain analytic and non-analytic pieces. 
The latter ones are genuine consequences of the spontaneous breaking of chiral symmetry, 
as they arise from quantum fluctuations or loops of the corresponding pseudo-Goldstone 
bosons (pions in the SU(2)p QCD case). They provide the complex contributions required 
by unitarity and come accompanied by UV divergences. Renormalization is ensured by the 
coefficients of the local terms in the expansion, the LECs. The values of these parameters are 
not determined by chiral symmetry and they enclose information on the dynamical content 
of the underlying theory, i.e. QCD. A consistent chiral power counting including baryons can 



be established considering the latter as non-relativistic fields j45|; for processes with a single 
baryon, a Feynman diagram with L loops, Nm meson propagators, A^ baryon propagators 
and Vk vertices of /cth order Lagrangian scales as 0(p n ) where 

n = AL - 2N M -N B + Y kV k- ( 4 ) 

k 

The chiral Lagrangian in the nucleon sector up to 0(p 3 ) can be expressed then as 12, |24| 

£ Ch PT = 4 2 2 + 4 4 2 + 41 + 4 2 i + 4 3 i , (5) 

where the superscript is the chiral order n. Here, C^J refers to the purely mesonic Lagrangian 
without baryons and 4lv corresponds to bilinears in the baryon fields. The explicit form 
for the mesonic Lagrangian is 

42 = ^>x + x + >, 

4 4 2 = ^4 (2<«X><X+> + <X+> 2 ) + • • • , (6) 
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where the ellipsis indicates terms that are not needed in the calculations given here and (■ • •) 
denote the trace of the resulting 2x2 matrix in the flavor space. For the different symbols, 
/ is the pion weak decay constant in the chiral limit and 



u 2 = U , u ll = iu'd^JJ u 



X± 



u'x u =t ux^u. 



(7) 



The explicit chiral symmetry breaking due to the non- vanishing quark masses (in the isospin 
limit m u = = rh) is introduced through x = IB^m. The constant B Q is proportional 
to the quark condensate in the chiral limit (m u = = 0), (0|g J V|0) = —Bof 2 8^. In the 
following we employ the so-called sigma-parameterization where 



U{x) 



ir{x\ 



7T (X) ■ T 



f 



For the pion-nucleon Lagrangian we have [24 



.(i) 

'ttN 
•(2) 



= Cl ( X+ )W - £^(u,u u )^D»D^ + h.c.) + |(uX># - jftf-fK, u u ]i/; + ... 



5m 



+- T^3 (K, [A,, n A ]](/W/J A + sym.) + h.c.) + ^([X-, + h.c.) 

C?14 — G?1 



,[£>„, u A ])£> A + h.c.) 



(9) 



In the previous equation m is the nucleon mass in the chiral limit and the covariant derivative, 
D^, acting on the baryon fields is given by = <9 M + T M with T M = [u\ o> M u]/2. On the 
other hand, the LECs (q and dj) are determined by fitting them to nN scattering data. 
Again, only the terms needed for the present study are shown in Eq. (Q. For further details 



on the definition and derivation of the different monomials we refer to Refs. 24, 46 



2.3. The Extended- On-Mass-Shell renormalization scheme 

The chiral power counting in the presence of baryons, Eq. (jl]), is blurred due to the 
presence of a new large scale, the nucleon mass m^, that is finite (does not vanish) in the 
chiral limit j2j and breaks the homogeneity of the amplitudes in the small scale p/A x . More 
precisely, in a Lorentz covariant framework, the loop contributions produce divergent pieces 
that violate the power counting formula requiring the renormalization of the lower-order 
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LECs. Any of these pieces will be denoted in the following as a power-counting breaking 
term (PCBT). 

The non-relativist ic treatment of the baryon fields leading to Eq. fll]) is implemented from 



the outset and in a systematic fashion within the HBChPT formalism [23j . As in the theory 
for mesons, the renormalization of the LECs in HBChPT can be completed order by order 
according to Eq. 01]), although at the price of loosing manifest Lorentz covariance. Moreover, 
as it was explained above, the non-relativistic expansion of the baryon propagators performed 
in the HBChPT scheme alters the analytic properties that a theory with dynamical nucleons 
should have, causing problems of convergence in some parts of the low-energy region. 

These problems can be overcome using modern manifestly covariant approaches. A 
consistent organization of the chiral power counting in covariant BChPT arises from the 
following crucial observation: The leading non-analytical behavior of the baryonic loop 



graphs obeys the power- counting formula (J3J) and agrees with the one given by HB [26|, |41 
This means that all the terms breaking the counting are analytical in quark masses and 
momenta and, consequently, of the same type as those given at tree-level by the most 
general chiral Lagrangian. A corollary is that one can trade the power-counting problem of 
the covariant approach for a renormalization prescription issue. 

The EOMS approach is a dimensional regularization scheme in which the bare LECs are 
adjusted to cancel the PCBTs present in the loop contributions. Notice that in this scheme 
one preserves exactly the right § analytic properties of the theory, in contrast with the also 
covariant IR renormalization scheme, in which a resummation of only the recoil corrections 
of the positive-energy part (particle) of the propagator is performed. To apply the EOMS 
scheme it is necessary to calculate analytically the terms coming from the loop integrals 
that can generate PCBT in the full amplitude. The technique we use is explained in detail 
Appendix D, and the EOMS renormalization of m, g and the 0(p 2 ) LECs is shown in 



m 



Appendix E 



2.4- The A(1232) resonance 

The resonances are an important feature of the low-energy hadronic spectrum and strong- 
interaction phenomenology. They appear as poles in the complex plane of the scattering 
amplitude or, more generically, of the QCD correlators. At low energies, the contribution 
of these poles can be expanded in Taylor series oip/5 around p = 0, where 5 is the scale 
of the mass gap between the ground state and the mass of the resonance. The role of the 
resonances in ChPT can be understood by means of a chiral effective field theory which 



includes not only the nucleons but all the resonances as dynamical degrees of freedom [47 . 
The conventional chiral Lagrangian is recovered by integrating out these resonances one- 
by-one and storing the information of their effects into the LECs. In fact, one can set the 
so-called Resonance Saturation Hypothesis (RSH) stating that all the short-range nature of 



the QCD interactions at low-energies is mediated by the resonances |47J, |48j. On the other 



hand, in cases where the mass gap 8 is not large enough, the p/5 expansion has a small 



3 By right analytic properties we mean those derived from S-matrix theory and implemented automatically 
in a (Lorentz covariant) quantum field theory of dynamical pions and nucleons. 
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radius of convergence, ruining the behavior of the chiral series. In this case, integrating out 
the resonance fields is not justified, so they ought to remain as genuine dynamical degrees 
of freedom in the theory. A remarkable example of this problem is the effect of the A(1232) 
in the ttN system for which 5 ~ — ~ 300 MeV, that is well within the expected 
region of validity of the chiral expansion. 

Treating the A(1232) as an explicit degree of freedom in a ChPT setup introduces two 
types of difficulties. The first one is related to the appearance of the new scale 5 and, thus, to 
the power counting associated to the Feynman diagrams including resonance lines. A power- 
counting method that takes into account the fact that M n < 5 < A x is the 5-counting 43j. In 
this method, one employs the power-counting formula Eq. (J3J) for the resonant contributions 
but using the assignment 5 ~ p 1 ^ 2 . This means that the A propagators, which count as 
~ 1/5, receive a suppression of 0(p 1 ^ 2 ) with respect to the nucleon propagators. Although 
in the (^-counting the expansion parameter is S/A x ~ 0.3, the convergence of the chiral series 
is expected to be faster than in the theory without the A. Notice also that this counting is 
only valid for the description of energies and pion masses below the scale 5, so it is not well 
suited for LQCD extrapolations or for describing the resonance region. In the latter case, 



the 5-counting changes to take into account the prominence of the A-pole [43|, [39 



The second issue is related to the representation of the resonances and the construction 
of suitable chiral Lagrangians. The A(1232) is a spin-3/2 resonance that can be described 
in terms of an isospin multiplet of Rarita-Schwinger (RS) fields A^ = (A ++ , A + , A , A - )^, 
where /i is the Lorentz index. The free RS Lagrangian is 

£ 3/2 = A M {i*r>>d p - m A7 n A„, (10) 

where 7 Miy and r y^ up are the anti-symmetric combinations of Dirac matrices. A very well 
known difficulty of a quantum field theory of high-spin particles concerns the consistency 
problem, which is due to the fact that the fields we use to represent these particles contain 
more components than physical degrees of freedom. For instance, the RS field is an object 
with 16 components (4 of the spinor times 4 of the Lorentz index) of which only 8 (4 in case 
of a massless field) correspond to the physical spin-3/2 particle and antiparticle. The Euler- 
Lagrange equations derived from Eq. ffTO]) provide the necessary constraints to guarantee, 
in the free case, that the independent components of the RS field are those corresponding 



to the physical degrees of freedom [49 



Introducing interactions that do not fulfill the right constraints upon the unphysical 



components of the RS field can lead to well known pathologies (see e.g. Refs. |49|, |39| and 



references therein). In order to tackle this problem, we adopt the consistent interactions, 
which are invariant under the transformation A M (x) — > A M (x) +d fl e(x). A remarkable prop- 
erty of these interactions is that they fulfill the same constraints as the free Lagrangian 



so they only account for the dynamical effects of the physical degrees of freedom [49]. A 



systematic procedure can be set to construct chiral-invariant La gran gians describing inter- 



actions among pions, nucleons and A's that are also consistent [39j, [50j. For instance, the 



leading 7riVA chiral coupling is given by [49 



4va = n~T~~ (<9pA M ) T^ u Nd^ a + h.c, (li; 



where T a are the spin-3/2 — > spin- 1/2 projectors, which verify T a T tb = S ab - r a T b /3. 
The parameter Ha is the NA axial coupling, which is poorly determined but related with 
the 7riVA coupling through the off-diagonal Goldberger-Treiman relation [5lj]. Using the 



one 



A(1232) decay rate, correspon ding to its Breit-Wigner width, = 118(2) MeV [52 
obtains the value ]%a = 2.90(2) [391 ] . 

At the order we work, higher-order 7riVA couplings can also contribute. In particular 
the following 0(p 2 ) chiral Lagrangians can be constructed 51 



£na = — NT a u a ^Y pa D P ^ - -y NT a u« ul ^D»D p A a + h.c, (12) 
where uj a = (t q [(9 m , u u ])/2. These Lagrangians are consistent and the on-shell equivalent of 



those accompanying the LECs 63 and b$ in Ref. [40 . 

For the sake of completeness, we also display the propagator of the A field, which is the 
ci-dimensional inverse operator of Eq. (llOp . 



p 2 — m\ + ie 



^ j 1p.1v 



d-2 

\1pVv ~ IvVy) - J-j TT-^PfiPv 



(13) 



(d-l)m A v "^ ^ (d-l)ml 1 

Notice that this propagator does not include the width of the resonance, which should be 
incorporated consistently order by order in the calculation of the specific Green functions. 

2.5. Isospin breaking 

Isospin symmetry implies a degenerate mass in the nucleon isospin doublet. However, 
the mass difference between the u and d quarks introduces a difference of 2.5 MeV in the 
neutron mass respect to the proton one. This splitting is compensated by electromagnetic 
interactions that rises the proton mass in 1.2 MeV. Thus, the overall correction to the 



neutron mass over the proton one due to the isospin breaking effects is of about 1.3 MeV [53 
Regarding the pions, the situation is reversed. In that case, the electromagnetic interaction 
produces the dominant isospin breaking corrections, while the quark mass difference starts to 
contribute at order (m u — rrid) 2 This gives a total splitting of the charged pions respect 
to the neutral one of approximately 5 MeV. These observations leads to the question about 
the impact of isospin violation in strong interaction processes like the one we are considering 
here. In fact, 7riV scattering is an excellent test ground for studying this question. High- 
quality data opened the possibility of studying the isospin breaking corrections from the 
experimentally accessible reactions ir + p and ir~p scattering and the single charge exchange 
(SCX) reaction n~p — > n°n. A study of the isospin and electromagnetic corrections within 



the framework of ChPT has been performed in Ref. [54J. In that paper, it is concluded that 
the maximal effects of isospin breaking in the irN scattering amplitude are approximately 
of the order of ~ —0.7% for the S- waves and of ~ —4% for the the P- waves, when CM 
energies below ^ 1.11 GeV are considered. In the present paper we work in the isospin limit 
and we will assume possible isospin breaking effects of this size when assigning errors to the 
PWAs that we use in Sec. HI 
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3. Calculation of the Scattering Amplitude 



In the following we display the details of the calculation of the irN scattering amplitude 
used in this paper. As it has been explained in the introduction, chiral representations of 
this amplitude have been obtained before in many different approaches and up to different 
degrees of accuracy. We present a ChPT analysis of the 7rjV-scattering amplitude up to 
0(p 3 ) accuracy that includes two main improvements over previous work. In the first place, 
we use Lorentz covariant B^PT in the EOMS scheme and, secondly, we explicitly include 
the A in the 5-counting and filtering the unphysical components of the RS spinors via the 
consistent couplings. We show in the next sections that this approach solves many of the 
problems reported by previous works. More precisely, it presents a natural convergence of the 
chiral expansion from the subthreshold region up to energies well above threshold, providing 
a suitable framework to undertake a model-independent analysis of all the phenomenology 
related to irN elastic scattering. We first revisit the chiral expansions up to 0(p 3 ) of the 
nucleon mass m^, the nucleon wave-function renormalization and the axial coupling g^, 
which are essential ingredients for the respective calculation of the ttN scattering amplitude. 

3.1. The nucleon mass, wave-function and axial coupling 

The chiral expansions of the nucleon mass and axial coupling are ubiquitous in baryon 



ChPT literature 15J. Expressions for and gj, in the EOMS scheme and up to 0(p 



accuracy can be found in Ref. |42|, |55| and [56|, |57| respectively. We take advantage of the 



re-derivation of these results up to 0{p ) to illustrate the application of the EOMS scheme in 
Lorentz covariant BChPT. For the nucleon mass, we need to consider the modifications of the 
nucleon propagator due to the term C\{x+)'4>ip m £>^n an< ^ the self-energy diagram in Fig. [2J 
An explicit calculation using the covariant Lagrangian and dimensional regularization gives 

* 



Figure 2: One-loop self-energy of the nucleon at 0(p 3 



the following expression for m^, 

m - AciM 
3g 2 M 3 



where 



m N 



3g 2 m . -. , , ,, u 3g 2 mM 2 
-fj r (2X(m 2 + M 2 )) + -» 



64tt 2 / 2 



rM 
— log 




(K) 


L m 


V 77T / 



4a 1 



M 2 
4m 2 



327T 2 / 2 



arccos 



2m) 



A 



m 



i[ln4vr + r , (l) + ll k 

16tt 2 U-4 2 L V ; 



(14) 
(15) 

(16) 



and where we have chosen \x = m for the sake of simplicity. As it was anticipated in Sec. 12. 3\ 
the loop contribution is 0(p 3 ) according to the power-counting formula (j4j) but it contains 
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analytic and divergent pieces which count as 0(p ) and 0{p 2 ). The EOMS method recovers 
the right hierarchy among the different contributions to mjy by choosing a suitable (and 
systematic) renormalization scheme, 



3g 2 m 3 .-. 
m — y m' H — — {2a), 



2/2 
3g 2 m 



(2A) 



3g 2 m 

1287T 2 / 2 ' 



leading to the final expression 



(17) 



m N = m! - 4c[M 2 



3g 2 M 3 rM 
647T 2 / 2 L m 



log 



4W1 



M 2 
4m 2 



arccos 



(M 
\2m 



18) 



One can explicitly check that this formula fulfills the power-counting and recovers the 0{p 3 ) 
HB result in a 1/m expansion, 



m7V = m ' - 4 C ;m 2 - ^ + e>( 



32tt/ 2 



A 2 m' 



(19) 



On the other hand, the wave function renormalization can be obtained straightforwardly in 
covariant ChPT, 



J N 



1 - 0>A 



3g 2 M 2 



327r 2 / 2 m 2 (4m 2 - M 2 ) 



2M(M 2 - 3m 2 WAm 2 - M 2 arccos 

\2m J 



+ {M 2 - Am 2 ) ((2M 2 - 3m 2 ) log I — J -2m 



(20) 




Figure 3: Diagrams that contribute to the nucleon axial form factor up to 0(p 3 ). The wavy line with the 
cross at the end corresponds to an external axial- vector source. The diamond coupling introduces the LEG 
di 6 . 
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The set of diagrams in Fig [3] contribute to the chiral expansion of axial coupling g& up 
to 0(p 3 ). A direct calculation in dimensional regularization gives, 



9a 



3 + 4M^ die - ^ ( - 2) m 2 + - l) M 2 ) A 



^M 2 



lQir 2 f 2 m 2 



P 

(3g 2 + 2) M 3 - 8 (# 2 + 1) m 2 M 



M 2 



arccos 



16/% 2 
M 



2m 



(3^2 + 2 ) m 2 + ( (4g 2 + 2) m 2 - (3/ + 2) M 2 ) log 



(21) 



where we have used again \x = m to fix the renormalization scale. In the first line, we have 
explicitly shown the analytic and divergent pieces which can be canceled by redefining the 
bare LECs g and d\§ in the EOMS scheme, 

, 2# (g 2 - 2) m 2 - g 3 m 2 



P 



-A 



16/ 2 7T 2 



16 



he 



g{g 7 



2f 2 



-A. 



'22) 



This leads to a renormalized expression of qa verifying the power-counting formula (Tjl), as i t 
can be explicitly seen by recovering the HB result, modulo an analytic ~ M 2 piece 




in the non-relativistic limit, 



g A = g' + 4M 2 d[ 



16 



g'M 2 

16n 2 



3g' 2 + (Ag 1 '' 



2) ^ {^) 

' \ m J 



0{ 



M 3 , 
A 2 m' 



(23) 



3.2. The ttN scattering amplitude 

According to the power counting formula 01]), the ttN scattering amplitude includes the 
diagrams of Figs.|U|5]and|6]at orders 0(p), 0(p 2 ) and 0(p 3 ), respectively. The contributions 
of the tree-level diagrams to the scattering amplitudes D ± (s,t) and B ± (s,t) are shown in 



the Appendix B The calculation of the loop diagrams is done by means of the Passarino- 
Veltman decomposition of tensor integrals in terms of 1-, 2-, 3- and 4-points scalar integrals. 
The results in dimensional regularization, and prior to renormalization, are presented in 
the |Appendix C| This decomposition of the loop-results provides a simple method to 



obtain the ultraviolet divergences and PCBTs of any of the loops that contribute to the 
amplitude in terms of those of the few scalar integrals (see Appendix D ). This facilitates 
the renormalization of the final scattering amplitude and the application of the EOMS 
scheme to preserve, manifestly, the power counting. 

The renormalization of the scattering amplitude is subtle. Only when all the loop dia- 
grams, including the pion wave function renormalization 



1 



6 V 



M 2 
~P 



2e: + 



167T 



log 



M 2 



+ 0{M$) 



(24) 
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Figure 4: Diagrams that contribute at LO, 0(p), to the irN scattering amplitude. The crossed it-channel 
has also to be considered. These diagrams enter at 0(p 3 ) with wave- function renormalization contributions 
and corrections coming from the fact we use physical values f„, gA and mjv (see text for details). 



o( P 2 ) 



Figure 5: Diagrams that contribute at NLO to the ttN scattering amplitude. The vertices denoted by a box 
introduce the 0{p 2 ) LECs ci_ 4 . 



and the nucleoli wave renormalization of the 0(p) diagrams in Fig. HJ are properly included, 
one can set a redefinition of the available LECs that casts the amplitude into a finite function 
fulfilling the power counting. The necessary re-definitions of the bare LECs in the EOMS 
scheme can be found in the Appendix [231 A remarkable consequence of chiral symmetry and 
the consistency of this scheme (as in any other <i-regularization scheme), is that the redefi- 
nition of the bare LEC c\ required to make the irN scattering amplitude finite and verifying 
the power counting is the same as the one demanded by the calculation of mjv, in Eq. ()17[) . 
and the pion-nucleon sigma term (see Sec. I5.3p . An important check of the renormalization 
and, thus, of our calculation has been to confirm that our results are independent of the 
renormalization scale@ This feature of the conventional dimensional-regularization schemes 
contrasts with the IR scheme, which introduces a spurious (higher-order) dependence on 

fJL 
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Our calculation is finally given in terms of the physical quantities f n , and gA, instead 
of their chiral-limit values which are not very well known. This procedure implies some 
reshuffling of 0(p 3 ) pieces. In particular, new contributions have to be considered arising 
from the 0(p) diagrams in Fig. H] when expressing the renormalized m, g and / in terms 
of the physical values m^, gA and f n . For the former two, we use the EOMS expressions 
derived in Eqs. ([18]) and (J2U [22]) respectively. For f n we use the conventional expression 



2 



M 

l + 7 r 



log 



167T 2 \ jJ? 



M 



0{Mt). (25) 



It is important to stress that we strictly keep the 0(p 3 ) contributions stemming from these 
re-definitions (e.g. in g 2 and m 2 entering the Born terms), so we avoid introducing any 



4 It is important to point out that some divergences at 0(p 4 ) are generated by the relativistic structure 
in the amplitude D + . This requires the inclusion of suitable 0(p 4 ) LECs, although we do not consider them 
explicitly. We first check that the residual scale dependence is negligible in a very wide range of values of fj, 
and then we remove by hand these divergences together with their associated scale dependence. 
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(a) (b) (c) 

\ / s / \ 




Figure 6: Diagrams that contribute at NNLO to the irN scattering amplitude. The vertices in the loop 
diagrams are 0(p) couplings, whereas in the tree level diagrams the diamonds introduce the 0(p 3 ) LECs 
d's. The diagrams in the crossed it-channel has also to be considered with the appropriate topologies. 



15 



Figure 7: Explicit contributions of the A(1232) resonance to the scattering amplitude up to 0(p 3 ) in the 
(5-counting. The crossed diagram is also included. 



high-order renormalization scale dependence. As a result, all the dependence on the LECs 
l\ and d[ 6 disappear from the chiral amplitude. On the other hand, the physical values can 
be directly used for the 0{p 2 ) and 0(p 3 ) scattering amplitude because chiral corrections are 
of higher order. 

For the numerical evaluation of the scattering amplitude, we programmed the scalar 



loop integrals and checked thoroughly the numerical results provided by LoopTools [58j . All 
the final numerical results and fits to PWAs are done with these subroutines. We use the 
numerical values / w = 92.4 MeV, M w = 139 MeV, m N = 939 MeV and g A = 1.267. 

Explicit contributions of the A(1232) appear up to 0(p 3 ) only through the Born-term 
diagram, Fig. [TJ and its crossed topology. At 0(p 3 ^ 2 ), one needs to include the terms 
given by the monomials in Eq. (ITT]) . At 0(p 5 ^ 2 ), Lagrangians in Eq. (fT2"|) have also to be 
considered. However, it has been recently shown that the latter couplings are redundant 
in the non-relativistic limit as their contributions can be absorbed in a redefinition of the 



coupling Ha and the LECs ci_4 59] . Nevertheless, we evaluate explicitly their contribution in 



order to numerically check the reach of these conclusions in a Lorentz-covariant framework. 



The results of these diagrams can be found in the Appendix B Finally notice that loop 
diagrams involving A lines start at 0(p 7 ^ 2 ) and so they are beyond the accuracy considered 
in this paper. 

4. Description of partial wave phase shifts 

In order to fix the LECs that appear in the Lagrangian and extract the physical informa- 
tion contained in the ttN scattering amplitude, we consider the ttN phase shifts provided by 



three different PWAs: The PWA of the Karlsruhe-Helsinki group |44j (KA85), the current 
solution of the George Washington University group fl (WI08), and the low energy phase 
shift analysis of the Matsinos' group j^] (EM06). 

The aim of partial wave analysis is to determine the hadronic phase shiftd^l from differ- 
ential cross sections and polarization data. They use the available database of 7r + p and ir~p 
scattering together with the SCX reaction n~p — > 7T°n. Assuming isospin invariance, all 
these reactions can be described by four invariant amplitudes, e.g. A ± (s,t) and B ± (s,t). 



5 In order to obtain the hadronic amplitude is necessary to take into account the Coulomb and electromag- 
netic corrections, which is are usually treated following [(30j], although EM06 implements these corrections 



in a different fashion (see Refs. pi IHlj). 
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One needs further theoretical constraints, in addition to unitarity (except near threshold), 
to provide a unique representation of the amplitude from the data [if. Such constraints 
are provided by fixed-t analyticity that, together with isospin invariance, are strong enough 
to resolve the ambiguities of the PWAs. On the other hand, there are still uncertainties 
resulting from experimental errors and discrepancies among different data sets (which are 
frequently more important). 

In what concerns the PWAs used here, both KA85 and WI08 assume fixed-t analyticity. 



While KA85 uses an old data set |62|, WI08 employs the set in [63|, which contains data 
collected in modern experiments. A very different approach is followed by EM06 which con- 
sists of a low-energy (-y/i < 1.16 GeV) phase shift analysis that employs hadronic potentials 
corrected electromagnetically jij]. In contrast with the methodology followed in KA85 and 
WI08, the Matsinos' analysis exclusively considers data in the elastic regionjfl Both EM06 
and WI08 have in common the inclusion of the new data collected along the last 20 years 
in meson factories. Studying these three different solutions provides a handle on systematic 
discrepancies that could arise from the different methodologies or the data sets employed in 
each of these analyses. 

After this brief introduction on the PWAs considered in this paper we proceed with the 



explanation of our fitting methodology. Our fits use Eq. ( Appendix A.ll ) for the calculation 
of the phase shifts, given that our calculation is perturbative. We also use the following x 2 

en(Si) 2 ' 

where 5 t h corresponds to the theoretical phase shift, while S and err(<5) denote the phase 
shifts provided by the PWAs and their errors, respectively. Since KA85 and WI08 do not 
give error for their phase shifts, we take for err(£) a sum in quadrature of a systematic error 



(e s ) and a relative error (e r ), as was done in Ref. [33(. For the systematic error we take 
the value e s = 0.1 degree, which is a value typically smaller than the value of the phase 
shifts. The introduction of a systematic error is advisable because avoids giving an excessive 
weight to the threshold region, where there is no experimental data. On the other hand, 
we take e r = 2% to take into account both the isospin breaking effects and the theoretical 
error coming from higher order corrections of 0(p 4 ), which are suppressed respect to the 
leading order by a factor (p/A x ) 3 . It is important to stress that a reasonable variation of 
these values does not affect the conclusions presented in this work. The specific choice we 
take is quite conservative if we compare with the errors assigned in the EM06 analysis. We 
consider two different approaches, the one without explicit A degrees of freedom (that we 
will denote as 4&-ChPT) and the one where we include the A as specified above (denoted as 
A-ChPT). 



6 It is worth mentioning that the Matsinos' group recently updated their PWA [64j |. The new results are 
very similar to the ones obtained in Ref. which are those employed in this work. 
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Figure 8: Fits to KA85 [44|] with $.-ChPT. The (red) solid lines correspond to the EOMS result 



and the (green) dashed ones to IR. Both fits are performed up to y/s n 



1.13 GeV. 



11. ^-ChPT 

First, we compare 4& results that stem from our calculation in the EOMS scheme with 
those obtained in HB Q and IR ^ at 0(p 3 ). We fit the phase shifts up to energies 
\J~s max = 1-13 GeV taking only the pion and the nucleon as the relevant degrees of freedom. 
In Figs. El M and [TU] we see qualitatively a good agreement between the EOMS and IR 
results for all the partial waves. On the other hand, in case of the fits to the EM06 solution 
(Fig. ITU]) , there is not published IR result to compare with. 

In Tabled] we show the values of the LECs resulting from the fits using the EOMS scheme 
(columns 2-4) together with the ones obtained using the IR prescription [33[ (columns 
5-6) and within the HBChPT formalism (last column). The fits are performed including 
the PWAs phase shifts up to y/s max = 1.13 GeV and the error quoted for the LECs in 
the EOMS results is determined, as in Ref. jUJ, by adding in quadrature the statistical 
uncertainty resulting from this fit to the spread of values obtained by varying the maximum 
energy ^fs max = (1.11, 1.12, 1.13) GeV in the fit. Notice that the fit to the EM06 solution 
has a very poor quality in terms of the x 2 ? and we only list the results at \fs max = 1-13 
GeV with the corresponding statistical uncertainty. This difficulty can be traced back to the 
problems of 4^-ChPT to describe the P33 partial wave (see Sec. 14.1. ip . The large value of 
the Xd o f i n this case affects considerably the theoretical prediction of the rest of the partial 
waves, in particular the P\\. 

From the comparison of the values listed in Table [TJ we can see that those determined in 
the EOMS scheme and using the different PWAs are compatible within errors. On the other 
hand, the values of the LECs extracted from a given PWA and using either the EOMS or 
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Figure 9: Fits to WI08 [4( with 4^-ChPT. The (red) solid lines correspond to the EOMS result and 
the (green) dashed ones to IR. Both fits are performed up to ^fs max = 1.13 GeV. 
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Figure 10: Fits to EM06 @] with jL-ChPT. Fits are performed up to y/s max = 1.13 GeV. 
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[331 
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3.72(32) 


3.66(31) 
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d\ + d 2 
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3.7(6) 


3.64(8) 


4.7(1.0) 


5.1(8) 


(3.2,4.1) 




-2.7(6) 


-2.6(6) 


-2.21(8) 


-3.9(9) 


-3.8(6) 


(-4.3,-2.6) 


d 5 


0.50(35) 


-0.07(16) 


-0.56(4) 


-0.3(5) 


-0.83(24) 


(-1.1,0.4) 


dn — ^15 


-6.1(1.2) 


-6.8(1.1) 


-6.49(2) 


-5.4(1.3) 


-6.2(1.1) 


(-5.1,-4.3) 


dig 


-3.0(1.6) 


-0.50(1.8) 


-1.07(22) 


-3.5(2.0) 


-1.1(1.5) 


(-1.6,-0.5) 


Xd.o.f. 


0.38 


0.23 


25.08 


0.45 


0.34 


(0.83 - 1.34) 



Tabic 1: Comparison between LECs in the different approaches of A-ChPT up to 0(p 3 ). The 
results of the present paper are given in the columns 2-4. The IR results 13311 are shown in columns 
5-6, and the last column corresponds to the ones obtained in HBChPT |24l ]. The errors shown for 
the EOMS case are obtained as it is explained in the text. 



IR are also consistent with each other and with those reported in HB 24(. Regarding the 
quality of the fits, in terms of the Xd f > ^ us stress that the EOMS and IR fits follow the 



same strategy, whereas the HBChPT result is taken from Ref. |24|, which uses a different 
error assignment. 

4-1.1. Convergence of the chiral series 

It is interesting to study the chiral expansion of the scattering amplitude calculated up to 
0(p 3 ) in the EOMS scheme by looking at the contribution of each order to the different PWA 
phase shifts. In Fig. [Til we plot the respective contributions to the total result (red line) for 
the A-ChPT fits to WI08. (Similar plots can be obtained for the KA85 and EM06 solutions.) 
This can be directly compared with Fig. 7 of Ref. (24| (HB) and Fig. 3 of Ref. [33[ (IR) 
since both are also A^ fits to PWAs in the isospin limit. As it is discussed in these references, 
we see that there exists a cancellation between 0(p 2 ) and 0(p 3 ) contributions in almost 
all the partial waves. Furthermore, the size of the 0(p 3 ) contributions can be very large 
and comparable to those given by the lower-order terms even at very low energies above 
threshold. Thus, the applicability of 4&-CI1PT to describe the PWA phase shifts at 0(p 3 ) 



is questionable. Nevertheless, the HBChPT study of Ref. 25| at 0(p A ) obtains that the 
corrections to the 0(p 3 ) result are more modest than those from the 0(p 2 ) to the 0(p 3 ) 
calculation, suggesting convergence. 

4.2. A-ChPT 

As it has been argued in Sec. I2.4[ integrating out the A(1232)-resonance from the chiral 
effective field theory is not well justified for the description of elastic irN scattering above 
threshold. In this section, we study this in more detail by applying A-ChPT within the 
EOMS scheme to the description of the PWA phase shifts. In the following, we present the 
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Figure 11: Convergence of the chiral series for the fit to WI08 in ^-ChPT. The short-dashed (blue), 
dashed (green), dotted (pink) and solid (red) lines correspond to the contributions of the 0(p), 
0(p 2 ), 0(p 3 ) and the total sum, respectively. 



results neglecting the 0(p 5 ^ 2 ) contributions of the 7riVA couplings d% and df (see Sec. 12. 4p . 
A discussion of the effects of these couplings is presented below. In Figs. [121 EH an d HH we 
show the results of the fits of A-ChPT to the KA85, WI08 and EM06 PWAs, respectively, 
following the same procedure for the fits as done in the theory. For the sake of comparison, 
we also plot the results previously obtained in 4&-CI1PT. In Table [2j we list the resulting 
values of the LECs. The central values stem from taking an average of the values obtained 
varying \fs max from 1.14 to 1.20 GeV (up to its maximum 1.16 GeV in case of the EM06 
solution), in intervals of 10 MeV. For the 0(p 3 ) LECs d's, we take averages weighted by the 
corresponding statistical uncertainty to reflect the fact that these parameters are accurately 
determined only in the fits at the higher energies (all the values from the fits at lower-energies 
are perfectly consistent within errors). The errors are obtained adding in quadrature the 
statistical error at \fs max = 1.20 GeV and the one resulting from the spread of central values 
at the different \fs maxl which, in general, gives the larger contribution to the error. 

As it can be deduced from these figures and the table, the inclusion of the A(1232) 
resonance has an important effect in the description of elastic irN scattering, given that 
now we are able to describe accurately the S- and P-wave phase shifts up to energies of 
v/i = 1.20 GeV and with xL.f. < L For the fits to the KA85 and WI08 PWAs, the 
Xdoi remains well below 1, whereas for the EM06 the improvement is more drastic since 
it reduces drastically the Xdof > f rom ~ 25 (4&-ChPT) to ~ 0.1 (A-ChPT). This illustrates 
the relevance of the A(1232) resonance even at very low energies and close to the threshold 
region. 
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Figure 12: Fits to the KA85 PWA[4Jj with A-ChPT (red solid lines) compared with the 4^-ChPT 



result (green dashed lines). The A-ChPT fits are performed up to \fs 



1.20 GeV. 
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Figure 13: Fits to the WI08 PWA[1] with A-ChPT (red solid lines) compared with the 4i-ChPT 
result (green dashed lines). The A-ChPT fits are performed up to y/s max = 1.20 GeV. 
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Figure 14: Fits to the EM06 PWAji] with A-ChPT (red solid lines) compared with the ^-ChPT 
result (green dashed lines). The A-ChPT fits are performed up to \fs max = 1.16 GeV. 
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c 4 
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di + d 2 
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-0.36(6) 


d 3 


-0.21(26) 


-0.23(27) 


0.28(4) 


d 5 


0.82(14) 


0.47(7) 


0.20(3) 


du — d±5 


-0.11(44) 


-0.5(5) 


0.35(9) 
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-1.53(27) 


-0.2(8) 


-0.53(12) 


h A 


3.02(4) 


2.87(4) 


2.99(2) 


Xd.o.i. 


0.77 


0.24 


0.11 



Table 2: Values of the LECs in A-ChPT. The errors are obtained adding in quadrature the statistical 
uncertainties at ^fs max = 1.20 GeV and the spread of values produced using different values of 
y/s max , from 1.14 to 1.20 GeV, in intervals of 10 MeV. The Xdof ^ s obtained at the maximum 
energies considered, yfs max = 1.20 GeV. 
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KA85 


WI08 

V V 1UU 


EM06 


RSH 


^1 


-0.46 


-0.50 


-0.47 


rv 


C2 


2.96 


2.73 


3.06 


1.9... 3.8 


c 3 


-3.78 


-3.59 


-3.91 


-3.8... -3.0 


C 4 


1.74 


1.65 


1.79 


1.4. . .2.0 


d\ + c? 2 


3.45 


3.55 


4.98 






-2.49 


-2.37 


-2.49 






-0.32 


-0.54 


-0.76 




C?14 — C?i5 


-5.99 


-6.30 


-6.84 




C?18 


-1.47 


0.30 


-0.54 





Table 3: Estimation of the A(1232) contribution to the 0(p 2 ) and 0(p 3 ) LECs in units of GeV -1 
and GeV~ 2 , respectively, based on the results obtained with (A-ChPT, Table [2]) and without {/^.- 
ChPT, Table [1]) the inclusion of this resonance as an explicit degree of freedom. In the last row 
we show the results of Refs. [I^] and 34] employing RSH. 



The effect of the A in the LECs is clearly shown in Table El In general, one notices an 
important reduction, in absolute value, of these parameters. As we will see below, these 
new values are more natural than the former ones as they lead to a good convergence of 
the chiral amplitude. It is interesting to compare the contribution of the A (1232) to the 
different LECs obtained explicitly from our fits with the one calculated using RSH 48| . An 
example of this comparison is given in Table [3] for the WI08 fits. We observe that there 
is a good agreement between our results and the RSH approach, except for c\, where we 
obtain that its value is shifted by an amount of ~ 0.5 GeV^ 1 . This can be interpreted as 
an indication of the fact that the LECs are stabilized once the tree-level A(1232) exchange 
contributions are taken into account (see also Sec. 14.2. ip . 

On the other hand, the results obtained fitting the different PWAs are grossly consis- 
tent with each other. Interesting differences can be found though, and these translate into 
discrepancies in the ttN phenomenology derived from the various analyses. Note also the 
larger error obtained in some LECs for the KA85 analysis. This originates from the nu- 
merical instabilities in the S"-waves of this solution, visible in the Fig. H21 and which lead 
to an over-estimation of the uncertainties for some observables (e.g. scattering lengths) 
derived from this PWA. The relative stability on the values of the LECs in our A-ChPT 
calculation contrasts very much with the results reported in A-HBChPT 40j, in which very 
large differences among different PWAs were reported, hindering a clear discussion of some 
related phenomenology, e.g. the pion-nucleon sigma term. For the value of the A(1232) 
axial coupling, tiA, we find that the WI08 solution gives a value that is perfectly compatible 
with the one directly extracted from the A(1232) Breit-Wigner width, Ha = 2.90(2). This is 
the width one should compare with because we are reproducing the phase shifts where one 
is sensitive to the Breit-Wigner shape of the resonance in the physical s-axis. For KA85, one 
obtains a coupling that is slightly larger, what could be related to the overestimation of the 
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width resulting from this PWA |44j. The EM06 solution also leads to a larger Ha than the 
one extracted from the width, although, we do not expect the EOMS analysis to describe 
accurately this quantity as it is focused in describing the phase shifts below the A-resonance 
region. It follows that the WI08 solution is the only one that gives a A(1232) Breit-Wigner 
width compatible with the value quoted in the PDG. 

We have included the tree-level A-contributions generated by the 0(p 2 ) nNA couplings 
in Eq. ( 1X2]) . As it was explained in Sec. [31 these couplings have been found to be redundant in 
the non-relativistic expansion since as they can be accounted for by a redefinition of the LECs 



and C!_ 4 |59j. We expect this reorganization of the chiral expansion to be effective also in 
the Lorentz covariant case because the leading contributions to the corresponding diagrams 
are their HB approximations. Nevertheless, we have checked this explicitly including the 
aforementioned pieces in the fits. We have found that: (i) the two LECs are extremely 
correlated with each other and only one of them can be kept in order to obtain stable fits. 



(Similar conclusions were derived in Refs. 24j and 5jJ.) (ii) Large correlations are found 



between the remaining df* and the LECs Ka and Cx_4. Therefore, the inclusion of the new 
LEC does not change appreciably the quality of the fits. (Hi) The values obtained for the 
df and different \/ s _ are stable and consistent with zero. In summary, we ratify the 
conclusions of Ref. |59|, what suggests that the chiral expansion can be organized also in 
covariant A-ChPT so as to remove these higher-order iVA couplings. Consequently, we do 
not include these contributions in the rest of the present work. 

Finally, we also studied the results of A-ChPT within the IR scheme. Although the 
EOMS and IR representations of the amplitude give a completely equivalent description of 
the phase shifts near threshold in the A case, they give very different results when trying to 
describe them up to the A(1232) region in A-ChPT. For instance, equivalent fits to those 
performed in the EOMS scheme at \fs max = 1-20 GeV lead to a Xd.o.f. °f 4.15 and 1.67 
for the KA85 and WI08 PWAs, that are much larger than the ones obtained in EOMS, 
0.77 and 0.24 respectively. Furthermore, the description of main observables related to 
the scattering amplitude, like the pion-nucleon sigma term, the GT discrepancy or the 
subthreshold coefficients, is not compatible with the results based on dispersive analyses 
which, on the other hand, are perfectly consistent with those obtained in the EOMS scheme 
(as we will see in the next sections). These difficulties arise from the IR- regularized loop 
contributions, which seem to develop some sensitivity to the u-channel unphysical cut at 
the energies reached in the A-theory (see also the discussion in Sec. 14. 3p . For completeness 



and in order to justify these conclusions we present in Appendix F a brief summary of main 
numerical results obtained in the A-ChPT within the IR scheme. Given these problems of 
IRChPT to describe the higher energies considered in A-ChPT we focus in the following on 
the representation of the nN scattering amplitude obtained in the EOMS scheme. 

4-2.1. Convergence of the chiral series 

It has been repeatedly argued in this paper that the chiral representation of the ttN scat- 
tering amplitude in A-ChPT presents better convergence than in the A-less case. Although 
the dramatic improvement in the description of the different PWA phase shifts achieved 
in the former approach indicates that this is the case, a conclusion in this regard can be 
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Figure 15: Convergence of the chiral series for the fit to WI08 in A-ChPT. The short-dashed (blue), 
dashed (green), dotted (pink), dash-double-dotted (orange) and solid (red) lines correspond to the 
contributions of the 0(p), 0(p 2 ), 0(p 3 ), A(1232) and the total sum, respectively. 



reached only after studying the contributions of the different orders to the amplitude. In 
Fig. [15l we show this comparison for the WI08 phase shift. Similar plots can be obtained 
for the KA85 and EM06 solutions. 

First of all, one confirms the expectation that the A(1232) is the main responsible for 
the rapid raise of the P33 phase shifts. Besides, the A-exchange gives a non-negligible 
contribution to the rest of the P- waves. Secondly, the contributions of the 0(p) pieces are 
significantly more important than the 0(p 2 ) ones in most part of the low-energy region and 
most of the partial waves (not for the Pn and P13 waves, for which these terms largely 
cancel each other). Nonetheless, the most important observation is that, in A-ChPT, the 
0(p 3 ) contributions are completely subleading compared to the LO and NLO terms in all 
the low-energy region above threshold. In terms of LECs, this means that their values in 
A-ChPT are smaller and more natural than in the 4&-CI1PT. Furthermore, in Sec. Owe will 
see that the A is also an essential ingredient to connect, in CM energy y/s, the subthreshold 
and the threshold regions. This whole picture in A-ChPT is consistent with the chiral power 
counting and it corresponds to what is expected from a well-behaved chiral expansion of the 
irN scattering amplitude. A remarkable consequence of these conclusions is that, while in 
4&-ChPT t he p oor convergence of the chiral expansion forces working in the subthreshold 



region [3J, |65|, A-ChPT can be applied to study the different phenomenology associated 
with the 7riV-scattering using, exclusively, the experimental information accessible in the 
physical region. 
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4-3. Unitarized amplitudes 

Another important comparison between the covariant schemes, EOMS and IR, stems 
from the resulting scattering amplitudes after taking the perturbative ones as input for non- 
perturbative S-matrix techniques. Proper unitarization approaches take care of the analytic 
properties associated with the right-hand cut and have proved to be very successful in the 
description of non-perturbative phenomena in ChPT 



66 



67] . In Ref. [33|, the IR representa- 



tion of the 7rN scattering amplitude was unitarized using an approximate algebraic solution 



to the N/D method that is obtained by treating crossed channel dynamics [66|, |67J, |68|. Al- 
though the covariant amplitude achieved a good description of the phase shifts, the unphys- 
ical cut introduced by the IR-method spoiled the description for energies yfs > 1.26 GeV. 
Given that the EOMS scheme has the right analytical properties, it is interesting to see if 
we can improve the description of the data, and explore the potential of the unitarization 
techniques applied on a reliable BChPT kerneL 

in order to compare with the IR 
in addition, a Castillejo-Dalitz- 



67. 68. 69 



We use the same unitarization method of 
approach Ref. [§3j|. For the P 33 partial wave we include, 

Dyson pole (CDD) [68| to take into account the contribution of the A(1232) resonance 
when unitarizing the amplitude of 4^-ChPT. In this method the unitarized amplitude, Tug, 
is written as 

1 
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s P 



for / ^ 3/2 or J ^ 3/2, 



for / = 3/2 and J = 3/2, 
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where Tui is the interaction kernel, 7 and sp are the residue and pole position of the CDD 
pole, respectively. The function g(s) corresponds to the unitary pion-nucleon loop, 
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where the subtraction constant is fixed by requiring that the Pn unitarized partial wave 
keeps the nucleon pole at the same position as in the perturbative calculation, i. e. physical 
nucleon mass. This condition translates onto the loop function as, 



g(s = m%) = 0. (28) 

On the other hand, we follow Ref. j69|, as it was done in Ref. 33], to extract the 
interaction kernel by matching the chiral amplitude obtained in 4&-CI1PT with the chiral 
expansion of Eq. (126]) order by order. Namely, taking into account that g(s) = 0(p), we use 
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T (l) + T (2) + T (3) = T (l) + T (2) + T (3) _ g ^ ( T (l))2 ^ (2g) 

where the superscript (n) refers to the chiral order of the amplitudes. Matching order 
by order we obtain the following relations between the chiral amplitude (T(s)) and the 
expansion of the interaction kernel T(s) 

rW( 5 ) = t«( s ), 
r (2) (s) = r( 2 )( s ), 

r (3) ( S )=T( 3 '( S )+ 9 ( S )(TW( S )) 2 . (30) 

In Figs. [16] and [T7] we show the fits of the unitarized amplitudes to KA85 and WI08 
PWAs (red solid line), and the previous result obtained within the IR scheme (green dashed 
line). We do not consider here the EM06 analysis because we consider data up to energies 
considerably higher than its upper limit. The EOMS fits are performed up to energies of 
\/s max ~ 1.3 GeV and achieve a very good description of data up to energies of \fs ~ 
1.35 GeV. The description of the WI08 phase shifts is better than for the KA85 solution, 
as it is reflected by a lower Xd.o.f. ( see Table HJ). Moreover, thanks to the CDD pole, the P33 
partial wave is described almost perfectly up to y/s = 1.35 GeV for both PWAs. Regarding 
the LECs, we see that fitting the unitarized amplitudes results in a set of values, for this 
range of energies, that is between the ^x- and A-ChPT perturbative results. This is not 
surprising since it is well known that the unitarization method used here respects the chiral 
order. Notice that by unitarizing the A-ChPT amplitude the values for the LECs have 
a smaller size than those in Tabled! except for d±s. On the other hand, the values for the 
different counterterms are very similar between the fits to KA85 and WI08 PWAs indicating 
that the procedure is stable. 

From Figs. [16] and [T7] it is also easy to see how the unphysical cut of IR (green dashed 
lines) affects the description of the phase shifts, giving rise to a sharp rise of the amplitude 
at around \/s » 1.26 GeV. This problem is clearly absent in the results obtained using the 
EOMS kernel, which has the conventional analytic properties. Although this issue arises at 
energies well above the expected range of ChPT, it sets a clear warning for future applications 
of the IR scheme in combination with unitarization techniques. This problem could be quite 
severe in the S77(3)-fiavor sector, where the effects of the unphysical IR cuts are already 
perceived in perturbative calculations [29 



5. nN scattering phenomenology 

Once the LECs have been determined, ChPT is able to predict a full set of related 
observables. In this section we study the results obtained for the threshold parameters, 
the Goldberger-Treiman relation and the pion-nucleon sigma term, in this order. We focus 
on the results obtained in A-ChPT although we also list those in 4&-CI1PT for the sake of 
comparing the two approaches. The results reported here for <j nN are already published 
elsewhere [7ol |. 
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Figure 16: Unitarized fits performed up to \fs max = 1.3 GeV to the KA85 solution. Solid (red) 
line: EOMS. Dashed (green) line: IR [331 ] 
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Figure 17: Unitarized fits performed up to y/s max = 1-3 GeV to the WI08 solution. Solid (red) line: 
EOMS. Dashed (green) line: IR Q 
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LEC 


KA85 


WI08 




^-UChPT 


^-UChPT 


Cl 


-1.04(2) 


-1.11(2) 


C2 


2.48(3) 


2.54(3) 


C3 


-4.48(5) 


-4.78(4) 


c 4 


3.00(2) 


3.04(2) 


d\ + d 2 


0.40(2) 


0.51(2) 


^3 


0.13(1) 


0.07(1) 


d 5 


0.53(2) 


0.36(2) 


du — d\s 


-0.33(4) 


-0.46(3) 


dis 


-3.90(13) 


-3.23(12) 


7 


0.0093(9) 


0.0072(6) 


s P 


2.38(7) 


2.21(5) 


Xd.o.f. 


2.71 


1.25 



Table 4: Value of the 0(p 2 ) (in GeV^ 1 ) and 0(p 3 ) (in GeV~ 2 ) LECs, together with the CDD 
parameters (in GeV 2 ), obtained in our best fits to the data of KA85 |44j and WI08 0). The 
label 4^-UChPT means that we took the perturbative result without the A(1232) to obtain the 
interaction kernel according to the technique described in this section. 



5.1. Threshold parameters 

The scattering lengths and volumes extracted from the different partial waves, together 
with the scalar-isoscalar (a^ + ) and scalar-isovector (a^+) scattering lengths, are shown in 
Table for the ^-ChPT and A-ChPT cases0 These results can be compared to the values 
reported by the PWAs S3, a a, which are listed in Table [6j One can see that there 
is a good agreement between the values extracted from the fits and the results of their 
respective PWAs. The only exception is the P33 scattering volume in the 4&-CI1PT fits. 
This is due to the A(1232)-resonance and to the inappropriate description of its effects in 
the 4^-theory j33j. In fact, we observe that the explicit inclusion of the A(1232) in our 
calculations improves the description of ap 33 . It is also remarkable that in A-ChPT, one 
obtains an accurate description of the threshold region despite that the fits are performed 
up to energies significantly above threshold. In the following we will focus on the discussion 
of the results obtained in this case. 

In the last column of Table we list the values of the scattering lengths obtained from 
pionic-atom data 0|. We can see that the PWA of the George Washington group, WI08, 
is the one that presents the best agreement with these independent experimental results. 
Regarding the results quoted for the Matsinos' group j^], it is important to point out that 
the values a£ + = 0.22(12)10~ 2 M" 1 and a^ = 8.78(11)10" 2 M- 1 are obtained by the same 



collaboration in studies of pionic hydrogen [72|. In this reference, only the scattering lengths 
for ix~p — > n~p and n~p — > n°n are provided and we calculate ad" + and Oq + using isospin 
relations. The comparison of the scalar-isoscalar scattering length Oq + is specially interesting 



7 For the technical methods employed to determine these parameters from the amplitude see Ref. [71( 
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Partial 


KA85 


WT08 


hlV106 


KAoo 


W108 


hlV106 


Wave 




A-^-hP 1 


a 7 nunT 1 


A-bnr 1 


A-LhP 1 


A-bhP 1 


+ 


— U.o(o ) 


U.4(o ) 


n c ( a \ 
U.D(4 ) 


-1.1(1.0) 


n 1 o/oo\ 
-U.L2{66 ) 


U.zo^zU ) 


a 0+ 


9.2(10) 


8.4(10) 


7.7(4) 


8.8(5) 


8.33(44) 


7.70(8) 


a S 3 i 


-9.9(13) 


-8.0(12) 


-7.1(6) 


-10.0(1.1) 


-8.5(6) 


-7.47(22) 


a Su 


17.5(21) 


17.2(21) 


15.9(10) 


16.6(1.5) 


16.6(9) 


15.63(26) 


a Pai 


-4.0(7) 


-3.5(7) 


-3.7(2) 


-4.15(35) 


-3.89(35) 


-4.10(9) 


a Pll 


-7.7(18) 


-6.0(18) 


-7.2(3) 


-8.4(5) 


-7.5(1.0) 


-8.43(18) 


a P 33 


25.1(9) 


23.7(9) 


23.6(2) 


22.69(30) 


21.4(5) 


20.89(9) 


a Pl3 


-2.7(7) 


-2.3(6) 


-2.7(3) 


-3.00(32) 


-2.84(31) 


-3.09(8) 



Tabic 5: Summary of the extracted values of the threshold parameters in ^i-ChPT and A-ChPT 
fits. The scattering lengths and volumes are shown in units of 10~ 2 M~ and 10~ 2 M~ 3 respectively. 



Partial 

Wave 


KA85 
[44J 


WI08 


EM06 


7r-atoms^ 


a o+ 


-0.8 


-0.10(12) 


0.22(12) 


-0.1(1) 




9.2 


8.83(5) 


7.742(61) 


8.71(10) 


a S 3 i 


-10.0(4) 


-8.4 


-7.52(16) 


-8.81(18) 




17.5(3) 


17.1 


15.71(13) 


17.5(3) 


a P . A1 


-4.4(2) 


-3.8 


-4.176(80) 






-7.8(2) 


-5.8 


-7.99(16) 






21.4(2) 


19.4 


21.00(20) 




a Pl3 


-3.0(2) 


-2.3 


-3.159(67) 





Tabic 6: Results for the threshold parameters obtained by the different PWAs (columns 2-4) and 
from 7r-atoms data in Ref. 0] (last column) . As before, the scattering lengths and volumes are 
shown in units of W~ 2 M ~ 1 and 10~ 2 M~ 3 , respectively. 



because this quantity is related to the not-very-well known scalar structure of the nucleon. 
In this sense, the result obtained from 7r-atoms data is perfectly compatible with our deter- 
minations of this quantity based on the modern WI08 and EM06 analyses. Besides, the error 
in the value of the KA85 solution is overestimated due to the unphysical oscillations of 
its S-wave data. Keeping this in mind, we conclude that our determination based on this 
older PWA is only compatible with negative values of the Qq + . 



8 Since the PWAs are based on ir + p and ir~p scattering data, the values given here are obtained from 
Table 6 of Ref. where the 7r + p and ir~p scattering threshold parameters are given in their isospin limit 
corrected by the isospin breaking corrections 73 1. 
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5.1.1. Chiral expansion of the threshold parameters 

In 1966, Steven Weinberg used current algebra to predict the ttN scattering lengths 



74j . and later, he re-derived them in the celebrated first application of non-linear chiral 



Lagrangians [191 ]. 

a+_ = 0(M 2 ), (31) 



In this section we investigate the corrections to these observables in ChPT |75( and we pay 
special attention to the convergence of their chiral expansions. The threshold parameters 
aue can be decomposed order by order in the following fashion 

„ - „« , J 3 /2) , (2) (3) mov 

aue - a Ije + a ul + a Ije + a Ije , [66) 

where we use the physical values of gA, fn and to calculate the leading-order terms. 
The Born-term, which counts as 0(p), gives a contribution to Oq + proportional to 0(M 2 ). 
We include this contribution into the 0(p 2 ) piece to avoid confusions with the Eqs. 032p . 
The numerical results of these expansions for and A-ChPT, and taking the results of 
the WI08 solution as an example, are shown in Table [71 In this Table we see that the 
convergence of the chiral expansion of the isovector scattering length is very fast in both 
approaches. In case of the isoscalar scattering length, the convergence seems to be slow in 
the A-theory and completely broken in the case. Besides that, the scattering volumes 
also present a problematic expansion in 4i-ChPT, as 0(p), 0(p 2 ) and 0(p 3 ) contributions 
are typically of the same size. This problem is alleviated considerably in the A-ChPT where 
a clear hierarchy 0(p) > 0(p 2 ) > 0(p 3 ) among the absolute values of the chiral corrections 
to each scattering volume is found. The A is important in these observables, in particular 
in the P33 volume for which its contribution is as large as the LO one. For the rest of the 
partial waves, the A contributions are subleading effects, of the same order as the 0(p 3 ) 
corrections. These results are consistent with the conclusions obtained from the inspection 
of the chiral expansions of the scattering amplitude discussed in Sec. HI 

5.1.2. Lattice QCD results 

In this section, we investigate semi-quantitatively the few results that have been reported 
by the LQCD community on the ir ^-scattering amplitudes. In particular, we study the 
results obtained using the quenched approximation by Fukugita et. al. j^j and also a recent 
calculation of the scalar I = \ phase shift near threshold for a pion mass of M n ~ 400 MeV 



reported by the NPLQCD collaboration [10j. The latter is specially interesting because the 
value given is large and has the opposite sign compared to the experimental results. A 
natural question that follows is whether the chiral extrapolation can explain the sign flip for 
this observable and, therefore, for the corresponding scattering length. In the following, we 
will use our 4^-ChPT results (LECs) to study the extrapolation to unphysical pion masses. 
Note that in this case we can not use the A-theory because the 5-counting is based on a 
hierarchy 5 > M n , that would be broken in the extrapolation. We take into account the 
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WI08 4^-ChPT 
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8 A-ChPT 






C?(p ) 


/■r\( 'i\ 

0[p ) 


bum 


0(p) 


/r\l 
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bum 







-0.48 


0.91 


0.43 








-1.04 


0.93 


n 1 1 r 

-O.llo 


°0+ 


7.91 





0.47 


8.38 


7.91 








0.42 


8.33 


a S 3 i 


-8.85 


0.46 


0.44 


-7.95 


-8.85 





-0.1 


0.51 


-8.44 




14.89 


0.45 


1.84 


17.19 


14.89 





-0.1 


1.77 


16.56 


a Pai 


-5.37 


3.90 


-2.08 


-3.54 


-5.37 


0.89 


1.03 


-0.44 


-3.89 


a Pll 


-18.15 


23.51 


-11.32 


-5.97 


-18.15 


3.57 


11.85 


-4.72 


-7.45 


a Ps3 


9.78 


13.11 


0.84 


23.73 


9.78 


7.67 


6.11 


-2.18 


21.38 


a Pia 


-4.89 


5.09 


-2.54 


-2.34 


-4.89 


0.95 


1.69 


-0.59 


-2.84 



Tabic 7: From left to right, numerical results for the chiral expansion of the threshold parameters, 
depending on whether the A-resonance is included or not, and using the fits to the WI08 PWA 
shown in Tables [2] and [H respectively. They are shown in units of 10 -2 M" 1 for the scattering 
lengths and 10~ 2 M~ 3 for the scattering volumes. The leading order coincides in both cases because 
we use the physical values of qa, fn and m^. 



M n (MeV) 


M p (MeV) 


U (MeV) 


M p 




732(5) 
527(6) 


989(8) 
876(11) 


144(1) 
120(2) 


0.74 
0.60 


0.40 
0.35 



Tabic 8: Results of Ref. [9( for M p (second column) and f w (third column), for different values of M w 
(first column). The fourth and fifth columns show an estimation for the maximum an minimum 
values of the expansion parameter in the chiral series for the corresponding pion mass. 



running of qa, fn and mjy with M n , using the expressions that has been given in Sec. [3] (we 
assume di$ = 0). 

To our knowledge, the work of Fukugita et al. jif is the only LQCD calculation of the 
7riV scalar scattering lengths (cto + and Oq + ) reported so far. It is important to recall here 
that using heavier-than-physical pion masses, as those required to connect with this LQCD 
calculation, necessarily slows down the convergence of the chiral expansion. The expansion 
parameter for describing the results of Ref. 0] can be estimated by taking into account the 
corresponding results for the p-resonance mass (M p ) and f n , taking naively that the chiral 
expansion scale lies between M p and 47rf n . Table [S] shows that the convergence of the chiral 
series at the pion masses employed in the LQCD calculation is, of course, poorer than at its 
physical value, where the expansion parameter lies between 0.12 — 0.20. 

With all this in mind, we show in Figure IT8lthe extrapolation of the threshold parameters 
Oq + and ag + from the physical point up to M n = 800 MeV for the 4&-ChPT result. The 
uncertainties of the LECs considered only include the statistical ones. These propagate into 
the extrapolations giving rise to the bands shown in the figure. We see that the fits to 
KA85 and WI08 are compatible, within errors, with the LQCD results. On the other hand, 
the EM06 analysis seems to disagree with the LQCD points. Nevertheless, the important 
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Figure 18: Results obtained from chiral extrapolations of the 4& case up to M w = 800 MeV. The 
errors taken for the LECs are taken from Table [TJ The points at the physical pion mass are taken 
from their corresponding PWAs, while the rest are taken from Ref. [jj. On the other hand, the 
green square at the physical pion mass shows the Weinberg's predictions. 



outcome of this exercise is to point to the fact that the extrapolation of the scattering lengths 
is very sensitive to certain combinations of LECs. Thus, LQCD calculations of the scattering 
parameters close to the physical point could be very important in the future to provide tight 
constraints onto the LECs relevant for ttN scattering and its related phenomenology. 



In fact, such a calculation has been recently reported by the NPLQCD collaboration [10 



In this work, the value of the Su phase shift at a CM energy 5E n ^ = 15.3 ± 1.8 ± 3.2 MeV, 
where SE nN is the energy respect to the nN threshold, has been extracted. Curiously 
enough, at that energy and M n = 390 MeV, the NPLQCD collaboration obtains a value of 
5n = —26 ± 7 ± 6 degrees for this phase shift, which is large and has a different sign as 
compared with the experimental result. The interest is then to see if BChPT can explain 
a sign-flip in the chiral extrapolation of this observable. We proceed as in the previous 
section and we extrapolate, for the Su partial wave, the result of the fits to the different 
PWAs, obtained with the physical pion mass, up to M n = 500 MeV and for an energy above 
threshold of 5E n N = 15.3 MeV. The results are shown in Fig. [19] and, as we can see, ChPT 
could, in principle, explain a change of sign for the scattering parameters associated to the 
.Su partial wave. As it also happened with the scattering lengths in Fig.[18]our results for the 
LECs fixed by fitting EM06 data are not compatible with the LQCD results. Nonetheless, 
it is worth recalling that the £ fits to the EM06 solution have a very large Xd o f an d are n °t 
very reliable. A closer inspection to the chiral structure of the observable under study shows 
that the chiral extrapolation depends on the cancellation produced between the combination 
of LECs C2 + C3, which is negative and — 2ci, which is positive. This is exactly the same 
cancellation occurring at 0(p 2 ) in the iso-scalar scattering length a^ + . Therefore, a sign flip 
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Figure 19: Extrapolation of the Su phase shift at energy SE^n = 15.3 MeV (see text for definitions) 
of the 4^-ChPT results. The point at the physical pion mass corresponds to the value given by 
each PWA. 



as the one described in Ref . 10] , could provide pertinent information on the scalar structure 
of the nucleon. 

5.2. The Goldberg er-Treiman relation 



The Goldberger-Treiman relation [17] is one of the most important and earliest results 
of the ideas of chiral symmetry applied to the strong interactions. It unveils a connection 
between a purely hadronic quantity, the pion-nucleon coupling g vN , and gA, which describes 
the hadronic coupling of the axial part of the weak current to the nucleon, 

gAm N 

9nN = — f — (1 + A GT ). (34) 

JTT 

This relation can be obtained directly using PCAC and the fact that the pseudoscalar current 
serves as an interpolating field for a pion. The breakdown of the relation, Aqt, is of order 
0(M%) due to the fact that these two couplings are evaluated at different kinematical points, 
with g n N at t = M% and g^ at t = 0. The smallness of the GT discrepancy inferred from 



this power-counting argument is ratified experimentally by studies based on NN PWAs [76 
and pion-atoms data [6j, leading to Aqt = 1 — 3%. 

In fact, the chiral expansion of Aq T does not contain non-analytic pieces up to 0(M%) (iif . 
A calculation of g n N in BChPT up to 0(p 3 ) proceeds essentially as it is shown in Fig. [31 
where the wavy lines are replaced by the pseudoscalar current coupled through a pion pole 
to the nucleonic lines. The non-vanishing contributions to the GT-discrepancy at this order 
only arise from the operator accompanying the LEC dis and from the loop correction of 
the type (d) in Fig. [3j that we label as A^ OMS . The Goldberger-Treiman deviation can be 
written as 

a _ 2M n d lS .(3) ,ox\ 

IX GT _ _ (- A EOMS , (doj 

9a 

where 

(3) . AMI / , Mt.fM, 



^MS=i^(l + 4>og(^)) +C W (36) 
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KA85 
^-ChPT 


WI08 
^-ChPT 


EM06 
4i-ChPT 


KA85 
A-ChPT 


WI08 
A-ChPT 


EM06 
A-ChPT 


Act 


10(4)% 


2(4)% 


3.7(7)% 


5.1(8)% 


1.0(2.5)% 


2.0(4)% 


9-kN 


14.1(5) 


13.1(5) 


13.29(10) 


13.53(10) 


13.00(31) 


13.13(5) 



Table 9: Results for Agt an d g n N from the different fits considered in this paper. 





KA85 

M 


WI08 


EM06 

u 


AW scattering 
[76J 


Pionic atoms 

m 


Act 


4.5(7)% 


2.1(1)% 


0.2(1.0)% 


1% 


1.9(7)% 


9-kN 


13.46(9) 


13.15(1) 


12.90(12) 


~ 13.0 


13.12(9) 



Table 10: Results for Agt and g^N obtained from the PWAs, AA-scattering and pionic-atom data. 



This result is exactly the same as in the relativistic calculation of Gasser et al. 0. The main 
contribution to A^q MS , of about 0.4%, comes from the 0(M%) piece, which in other schemes 
such as IR and HBChPT, is absorbed through a redefinition of di$. The 0(M%) pieces 
from Eq. fl36|) . according to a naive power counting estimate, contribute only with a tiny bit 
~ 0.01%, suggesting that the chiral expansion in Aqt converges at very fast pace. Therefore, 
Agt is described, in very good approximation and with a negligible theoretical uncertainty, 



by Eq. f[3"5]) [251 ] . The pion-nucleon coupling g nN , and the respective GT-discrepancy, can 



be extracted from the ttN scattering amplitude using numerical methods related to the 



determination of the residue of the amplitude at the nucleon pole |33l. I71l|. We have checked 
that the application of these methods to the EOMS-renormalized amplitude naturally lead 
to the same values of g n N as those obtained by the direct application of Eq. (1351) . 

In Table [H we display the results on g n j^ and Aqt obtained from the different fits in 
the EOMS scheme considered in this work. The uncertainties are propagated from the 
errors in d 18 shown in Tables [TJ and [2J In Table [TD1 we show also the values reported by 
the different PWAs and the ones that are obtained from independent experimental sources, 
AA-scattering [76[ and pionic-atom data pj. First of all, notice the larger errors in the 
ChPT values, due to the troublesome convergence and, correspondingly, larger uncertainty 
in the determinations of the LECs in this approach. Also, it is worth pointing out the larger 
error in the A-ChPT value obtained for WI08 results due to the sensitivity to di$ in the fits 
to this PWA when varying \fs max - 

The comparison between both tables shows that the determinations from the fits of 
A-ChPT to the PWA phase-shifts are quite consistent with the values reported by the 
respective collaborations. Interestingly enough, we see that not all the PWAs are consistent 
with the extractions of g n ^ obtained from alternative experimental sources. In particular, 
we see that the WI08 and EM06 results agree with those obtained from AA and pion- 
atom data. On the other hand, the KA85 value is considerably larger than any of the 
other determinations, leading to a GT-discrepancy with a size that is currently considered 
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implausible [34J . Therefore, the analysis of g^N gives phenomenological support to the WI08 
and EM06 solutions, in detriment to the KA85's one. 



5.3. The pion-nucleon sigma term 

The definition of the pion-nucleon sigma term can be given in terms of the commuta- 



tors 



77 



= 2^T ^ I {N(p) 1 [Q% [Q% Hsb]] \ n{p)) ' (37) 

a=l 

with Q a A the axial charge and T-Lsb is the chiral-symmetry breaking part of the QCD Hamilto- 
nian. This quantity is related to the explicit breaking of the chiral symmetry and, therefore, 
it should be small compared to A x . From Eq. ( 137]) it is straightforward to obtain 

77? — 

= „ (N(p)\(uu + dd)\N(p)), (38) 

where m = (m u + m^)/2. The sigma term, in this form, can be identified with the nucleon 
scalar form factor u(p')a(t)u(p) = (N(p')\m(uu + dd)\N(p)) evaluated at t = 0. This matrix 
element can be derived also by means of the Hellmann-Feynman theorem from the quark- 
mass dependence of the nucleon mass, 

o dm ]v 

**» = M "q^- ( 39 ) 

An explicit calculation of the scalar form factor or the nucleon mass in EOMS up to 0(p 3 ) 
gives 

a nN = -4c M - -MA ■ * N ' arccos + M„ log -1 , 40 



where c' x is the LEC renormalized in the EOMS scheme (see Appendix E ). 

As it was discussed in the introduction, the pion nucleon sigma term is a quantity related 
to the structure of the nucleon which is important to understand the origin of the mass of 
the ordinary matter and the nature of the breaking of chiral symmetry in QCD. An accurate 
value of this matrix element is also required in to reduce the hadronic uncertainties that enter 
in the phenomenology of direct searches of dark matter. The main method to determine a n ^ 
experimentally is by analytical continuation of the isoscalar scattering amplitude to the CD 



point [18[ (see also Sec. 16. ip . However, there is still no consense on the value of oy/ v because 



it varies depending on the PWA taken as input [78|, [79|, |80] • In this respect, Ref. [70[ made an 
important step forward by obtaining perfectly compatible results between ChPT and those 
reported from the dispersive analyses based on the same PWA. It should be stressed that 



both methods are well sound and model independent. The agreement reported in Ref. [70 
makes then clear that the problem to fix cr^/v rests on the data basis employed and not on 
theory. 
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KA85 


WI08 
A-ChPT 


EM06 
A-ChPT 


KA85 
A-ChPT 


WI08 
A-ChPT 


EM06 
A-ChPT 


a wN (MeV) 


79(11) 


97(7) 


95(3) 


43(5) 


59(4) 


59(2) 



Tabic 11: Results for a^N in MeV for the cases without and with the A(1232), A-ChPT (Tabled]) and 
A-ChPT (Tabled)), respectively. 





KA85 [78 


WI08 [79 


EM06 [80; 


a nN (MeV) 


45(8) 


64(7) 


56(9) 



Table 12: Results for ova?, in MeV, extracted by the different PWAs. 



In this section we briefly review the results presented in [70j. The chiral Ward identity 
relating the scattering amplitude and a n N a t the CD point, can be accessed in a more 
elegant way using ChPT. Indeed, a value of a n N can be predicted once the relevant LEC c\ 
is properly determined from the scattering data. However, a reliable value for a n N would 
only follow from a representation of the scattering amplitude with a well behaved chiral 
expansion above threshold. 

In Table [TTj we list the results for obtained for the different PWAs. The results 
derived from the 4^-ChPT are shown for completeness. In Table [12] we also give the results 
obtained from dispersive analyses and using the CD theorem. As we can see, the results 
obtained in the A-theory are quite accurate and perfectly consistent with the dispersive 
results. On the other hand, the values obtained for the 4^-case tend to be larger. As 
it was discussed in Sec. 14.21 this sizable effect of the A in the values of c\ and a n N is 
not expected on the grounds of the RSH. Finally, we reported in Ref. 70[ the value g^n = 



59(7) MeV as it is extracted from the WI08 and EM06 PWAs. The error includes systematic 
uncertainties added in quadratures to the theoretical uncertainty that is estimated from the 



explicit calculation of higher-order diagrams. See Ref. [70j for a detailed analysis of this 
determination and for the discussion of the consequences that a relatively large value of 
ct^tv have in phenomenology. In relation with the strangeness content in the nucleon we 
address the reader to the recent reanalysis of the so-called "strangeness puzzle" presented 



in Ref. [81| . There it is shown that the relatively large value reported here for g^n is not at 



odds with a small strangeness content in the nucleon. 



6. Subthreshold region 



A proper description of the subthreshold region is very important in BChPT because 
the so-called soft point, s = u = m 2 N and t = 0, is the one about which the chiral expan- 
sion is performed (this point is defined in the chiral limit). Previous BChPT analyses have 
found difficulties to connect the information around the soft-point (as derived from disper- 



sive studies) with the experimental data [24j, |25|, |34j. As a result, it has been concluded 



that BChPT at one loop is not accurate enough to relate the subthreshold and physical 



38 







T/" A O C 


\T 7TAO 
W1U8 


TT11\ f fir 1 


J\Ao5 


W1U8 


tiMUb 






/V / 11 . r>rn 




A-ChP 1 


A-ChP 1 


A f ^ l . T> r n 

A-bnr 1 


A-bnr 1 




/" /\/f— 3\ 
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a 01 


i. <3(19 j 
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u 02 
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0.021(6) 
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0.005(6) 


0.004(6) 
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-5.1(1.7) 


-5.1(1.7) 


-4.5(9) 


^00 




1.81(24) 


1.68(16) 


1.495(28) 


1.63(9) 


1.53(8) 


1.379(8) 


d m 




-0.17(6) 


-0.20(5) 


-0.199(7) 


-0.112(25) 


-0.115(24) 


-0.0923(11) 


d w 


(M- 2 ) 


-0.35(10) 


-0.33(10) 


-0.267(14) 


-0.18(5) 


-0.16(5) 


-0.0892(41) 




17(7) 


17(7) 


16.8(7) 


9.63(30) 


9.755(42) 


8.67(8) 



Tabic 13: Results for different subthreshold coefficients obtained from the LECs shown in Tables [T] and [5] 
obtained from fits to the PWA phase shifts in fa- and A-ChPT, respectively. 



regions 34J. In this section we study the extrapolation of the chiral representation of the 
scattering amplitude into the subthreshold region in BChPT within the EOMS scheme. We 
also revisit important low-energy theorems established at specific points of the subthreshold 
region: The CD point 18] and the Adler point (2ol |. 

The starting point is the so-called subthreshold expansion, 



'10' 



"20 1 



c±t 2 + ... 



(41) 



with X ± 
that v = 



0. 



D- 
t -- 



D /v,B + /is, B , the Born-subtracted scattering amplitudes [l|. Notice 



corresponds to s 



u 



m 



N 



+ M 2 , so that the expansion in Eq. ( 14 ip 



is done around the same point as the chiral expansion but for physical pion masses and 
with the coefficients xfj non-analytic functions of M w . The usual procedure to analyze the 
subthreshold region of the ttN scattering amplitude is to determine the values of the leading 
coefficients in Eq. ( l4T|) using dispersive analyses [lj]. The resulting values can then be used 
to fix the LECs of BChPT and to analyze the consequences of chiral symmetry around the 



soft point (low-energy theorems) [65|, |34 



In this work we follow the inverse procedure as we fix the LECs in the physical region, 
where data actually exists, and we investigate the resulting description in the subthreshold 
region. In Table [13] we show the results for the leading subthreshold coefficients obtained 
from the fits to the PWA phase shifts in 4&-CI1PT and A-ChPT. In Table [T41 we list the results 



for the same coefficients as they have been reported by the different PWAs [44j, |79|. The 
first thing that is worth noticing from the comparison between the two tables is to confirm 
that z^-ChPT in the EOMS scheme fails to connect the physical and subthreshold regions 
at 0(p 3 ). Indeed, the numerical values of the subthreshold coefficients are not consistent 
with the dispersive results, even though the ChPT values have a sizable uncertainty. This 
problem is not likely to be solved at 0(p A ) in z^-ChPT in EOMS in the light of the results 



obtained in HBChPT 25 



and IR [34J at this order. 
On the other hand, A-ChPT gives a description of the subthreshold region that is, in 
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-0.134(5) 




d w 
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-0.167(5) 




b oo 


10.36(10) 





Table 14: Results on subthreshold coefficients from the Karlsruhe and George Washington groups. 



general, perfectly consistent with the dispersive results. The only disagreement concerns 
the coefficient g^}~ 2 , which has its physical origin in the incapacity of BChPT to reproduce 



properly, at this order, the curvature induced by the two-pion threshold at t = 4M 2 l78| 
(This has a very important consequence on the determination of a n N using the value of the 
scattering amplitude at the CD point, as we will see in the next section.) We conclude, 
then, that the explicit A-exchange contribution is a fundamental ingredient to bridge the 
gap between the physical and subthreshold regions. This is an important result as it paves 
the road for studying all the phenomenology related to irN scattering in a systematic manner 
within A-ChPT, using directly scattering data and without any other dispersive input. This 



contrasts with the the conclusions derived in A-HBChPT |40j. On the other hand, further 
studies at higher orders in the chiral expansion and including the A (1232) degrees of freedom 
in a coherent way should corroborate this particular finding of the present work. Such a self- 
consistent framework to study the ttN scattering amplitude, based exclusively on BChPT, is 
complementary to other model-independent approaches based on a pure dispersive treatment 



of the amplitude, e.g. [82 



Finally, it is also interesting to compare the results obtained for the subthreshold coeffi- 
cients from the different analyses and in A-ChPT. The KA85 and WI08 results closely agree 
with each other, besides an important discrepancy in o^ which is related to the different 
pion-nucleon sigma terms reported by the two solutions. However, comparing those with the 
novel results obtained for EM06 solution, we see that the latter gives a physical picture of 
the subthreshold region around the point {y — 0, t = 0) that is quite different to the former 
ones. In fact, the values of most of the subthreshold coefficients obtained from the fits to the 
EM06 phase shifts are not compatible with the ones extracted from the KA85 or the WI08 
solutions. The fact that the two latter PWAs grossly agree gives support to their solution 
in the subthreshold region. This discussion could also take place at the level of the values 
of the LECs (as shown in Table [2]), which the subthreshold coefficients ultimately depend 
on. Nevertheless, it is worth remarking that meaningful comparisons among different PWAs 
can only be done based on observable quantities. 
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KA85 
A-ChPT 


WI08 
^-ChPT 


EM06 
^-ChPT 


KA85 
A-ChPT 


WI08 
A-ChPT 


EM06 
A-ChPT 


% 


-16(18) 
1.0(1.1) 


4(11) 
0.3(8) 


6.0(1.7) 
0.40(11) 


-17(8) 
1.0(5) 


-4.4(4.5) 
0.28(28) 


7.3(2.1) 
0.50(14) 



Tabic 15: Check of the Adler condition. The second and third row show the deviation from this 
condition, in MeV, and the relative value of this deviation, respectively. 



6.1. The Adler consistency condition and the Cheng-Dashen theorem 

The isoscalar scattering amplitude D + (i>, t) is subject to a couple of important low-energy 
theorems. Its extrapolation onto t > at v = is constrained by the Adler consistency 



condition 20| at t = M% and relates the amplitude to the pion-nucleon sigma-term at the 
CD point, t = 2M% 18[. The Adler's consistency condition states that 



D + (u = 0,t = Ml)~^, (42) 

m N 

which is equivalent to the statement that the Born-subtracted isoscalar amplitude has a zero 
in the neighborhood of t = M%, D + {v — 0, t = M%) ~ 0. Once the LECs have been fixed by 
fitting our theoretical amplitude to the PWAs data, it is interesting to check the consistency 
of our BChPT calculations with the Adler's condition. 

In Table [T5l we display the results for f%D + (v = 0, t — M%) and the different fits per- 
formed in this work. We also list the relative deviation with respect to the exact fulfillment 
of the Adler condition, Eq. (|42]) . As shown in this table, there is no much difference between 
the results obtained for A-ChPT and A-ChPT, that is not surprising since the Adler con- 
dition is a direct consequence of PCAC and so it has to be satisfied in both cases. On the 
other hand, the WI08 and EM06 analyses fulfill the Adler condition better, although KA85 



is also within the theoretical bounds established by Adler [20(. Finally, these results can 
be compared with those given by the dispersive calculation, which can be obtained using 
Eq. (j4"Tj) and the respective values of the subthreshold coefficients, leading to —16 MeV and 
—4 MeV for KA85 and WI08 solutions, in order. These values compare well with the ones 
obtained in both versions of BChPT. 

The second low-energy theorem is more important as it relates the isoscalar amplitude 
to the pion-nucleon sigma term at the CD point, 

S = flD + {v = 0, t = 2Ml) = a wN + A CT + A R . (43) 

In this equation, D + (u,t) is the Born-subtracted isoscalar vrA^ scattering amplitude, A CT = 
it(2M|) — cr^N and Ar is a remainder originating from the translation of the exact relation 
at the soft-point (and thus with off-shell pions) to the CD point [l8[. The different pieces 
appearing at both sides of this equation can be obtained in BChPT 0, |83[. The difference 
of the scalar form factor at t — 2M% and t = in the EOMS scheme gives the same result as 
the one obtained by Gasser et al. in their seminal paper jij, A CT ~ 4.7 MeV. The remainder 
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KA85 
^-ChPT 


WI08 
^-ChPT 


EM06 
A-ChPT 


KA85 
A-ChPT 


WI08 
A-ChPT 


EM06 
A-ChPT 


E (MeV) 


84(11) 


103(7) 


103(3) 


48(5) 


64(4) 


64(2) 



Table 16: Results, in MeV, for the E-terms obtained in our analysis. 



of the CD theorem, Ar, is an analytic piece of order 0(M%) @, 0, which comes out to 
be numerically very small Ar < 1 MeV In Table fT6| we show the results for E in the 
different schemes treated in this work which can be compared with those of a n N in Table [TTJ 
in order to confirm the fulfillment of the CD theorem in our calculations. 

However, it is known that, at 0(p 3 ), BChPT fails to catch the full strength of the two- 
pion threshold in the extrapolation of the scalar form factor of the nucleon to t = 2M% [781] Fl 
Namely, dispersive calculations of the difference between these two kinematical points lead 



to the value A^ ~ 15 MeV [78|. This issue could limit the applicability of the CD theorem at 
0(p 3 ) as the determination of a n N from the 7riV scattering amplitude might be afflicted by 
a systematic uncertainty of ~10 MeV. On the other hand, comparing the A-ChPT results 
in Table [16] for the KA85 and WI08 solutions with those given by the dispersive analyses, 



E ~ 60 MeV |78J and E = 79(7) [79j MeV respectively, indicates that the same problem 



afflicts the t-dependence of the scattering amplitude on the left-hand side of Eq. f T43 



Namely, the quantity E can be rewritten as [78 



E = S d + A D , (44) 

where Ej = f^jdtn + ^M^d^) and Ad is the remainder given, in very good approximation 
(SH ~ 1 MeV [T9|] ) , by the curvature term, A D = 4f%M*d£ 2 . Neglecting A R , the CD theorem 
now takes the form 

a wN = E d + A D - A a . (45) 

In the previous section, we found that the values obtained in A-ChPT for the first two 
terms of this expansion (d^Q and dpi) agreed with the ones extracted from the dispersive 
analyses, so there is also agreement on the determination of E^. All the discrepancy between 
the values of E extracted in A-ChPT or the dispersive analyses thus originates from the 
discrepancy on the d^ 2 coefficient. As it was discussed above, this is related to the fact that 
the scattering amplitude at 0(p 3 ) does not catch the full strength of the two-pion threshold, 
which translates into an underestimation of ~ 10 MeV in the value of Ad and, hence, of E, 
when extrapolating to the CD point. 

In conclusion, the crucial point in the determination of a^N is not the value of E but 
rather the one of Ed, which is properly given by A-ChPT, together with the value of the 
difference Ad — A CT . In this sense, the very same effect curving the ^-dependence of the 
isoscalar scattering amplitude, enhances the slope of the scalar form factor at t — such 



9 This effect is enhanced by the particularly large scalar isoscalar irir partial wave amplitude driving the 
scalar form factor of the pion . 
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that these contributions largely cancel in the extraction of the sigma term. Indeed, although 
A-ChPT fails to give a reliable description of Ad and A CT individually, it gives a value for 
their difference, Ad — A CT = —3.5(2.0) MeV, which is perfectly consistent with the dispersive 
result A D — A CT : 



-3(1) MeV 78 



This analysis clarifies the reason why the pion-nucleon 
sigma term can be accurately pinned down in A-ChPT at 0(p 3 ) with the residual uncertainty 
produced by these effects well within the theoretical error estimated through the explicit 



calculation of 0(p 7 ^ 2 ) and 0(p 4 ) diagrams 70 



7. Summary and conclusions 

We have presented a novel analysis of the irN scattering amplitude in Lorentz covariant 
B%PT within the EOMS scheme up to 0(p 3 ) and considering the inclusion of the A(1232) 
explicitly in the ^-counting. We first studied the phase shifts in partial waves provided by the 
Karlsruhe-Helsinki, George- Washington and Matsinos' groups, that we use as experimental 
data to fit our LECs. While the ^-ChPT approach has the same difficulties to describe 
the region above threshold as those found previously in the HB and IR schemes, the A- 
ChPT perfectly describes the phase shifts up to energies below the A-resonance region. The 
improvement achieved in the latter case is clearly illustrated by the analysis of the EM06 
solution, for which the 4& fit gives a very large Xd o f > that ^ s P u ^ we ^ below 1 once the A is 
explicitly included. Differently to a previous analysis up to the same accuracy in HB and 
SSE scheme, the values of our LECs are stable against the PWAs phase-shifts used as input, 
allowing for a clear discussion of different irN phenomenology. 

Once the LECs are determined, we study thoroughly all the observables associated with 
the ttN scattering amplitude. In particular, we discuss the results and chiral expansion of 
the threshold coefficients, the Goldberger-Treiman relation, the pion-nucleon sigma term and 
the extrapolation onto the subthreshold region. Also, we investigated semi-quantitatively 
the extrapolation of the scarce results on scattering lengths and Sn phase shift reported by 
the LQCD community. In general, we conclude that the A-ChPT converges much better 
than the A approach and that, in the former case, one obtains a phenomenology perfectly 
consistent with the one reported by the PWAs. From the comparison among the results 
on observables that are obtained by an analysis of the different PWAs, we conclude that 
the WI08 solution is the most consistent with those extracted from alternative experimental 
sources (iViV-scattering and pion-atom data). We remind here that the KA85 analysis gives 
rise to a value for Ha that is not compatible with the value obtained from the A (1232) 
Breit-Wigner width (in agreement with the KA85 overestimation of this observable) and 
to a value for g n N that leads to a sizable violation of the GT relation, which is nowadays 
theoretically implausible. As for our study of the EM06 PWA, we found a value for the 
isovector scattering length that is too small as compared with the accurate values obtained 
from pion-atoms data. Besides, the picture of the subthreshold region arising from this 
solution around the point {i> — 0, t — 0) is quite different to the ones given by KA85 and 
WI08. 

The most important conclusion of our work is that the scattering amplitude in A-ChPT 
converges well from the subthreshold region up to energies well above threshold. This shows 
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that a systematic framework to analyze the irN elastic scattering data without spoiling the 
structure of the amplitude in the subthreshold region is possible using the A-ChPT approach 
developed in this paper. This is a remarkable result since it should allow to extract all 
the observables related to the irN scattering amplitude directly from the differential cross 
sections at low energies and in a completely model-independent fashion. Also, a calculation 
up to a 0(p 7 / 2 ) in A-ChPT is called for to confirm the good behavior of the chiral series in 
our approach and to improve the theoretical uncertainties of our determinations. 
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Appendix A. Partial Wave Decomposition 

The free one-particle states are normalized according to the Lorentz-invariant normal- 
ization, 

(p', a'; 7|p, a; 7) = 2£ p (2tt) 3 <5(p / - p)<W<5 7 y (Appendix A.l) 

where E p is the energy of the particle with three-momentum p, a the spin of the nucleon 
and 7 indicates any internal quantum number. A free two-particle state is normalized 
accordingly and it can be decomposed in states with well defined total spin S and total 
angular momentum J. In the CM frame one has, 

\n(~p;a)N( Pl a;a)} = V4^J2(mafi\£SJ)Y e m (py\Jfi£;aa), (Appendix A.2) 

with p the unit vector of the CM nucleon three-momentum p, a and a the isospin third- 
components in the Cartesian basis of the pion and nucleon, respectively, £ the orbital angular 
momentum, m its third component, fi — m + a the third-component of the total angular 
momentum and S the total spin, with S = 1/2 for irN scattering. 

The Clebsch-Gordan coefficients are denoted by (mim2m3|jij2j3), corresponding to the 
composition of the spins j\ and j'2 (with third-components rrti and m 2 , in order) to give 
the third spin j 3 , with third-component m 3 . The state with well-defined total angular 
momentum, \Jfi£;aa), satisfies the normalization condition, 

4:71 /S 

(J'/j,'£';a'a'\Jiil)aa) = SjjrSa'uSu'—r-rSa'aSa'a- (Appendix A.3) 

IpI 

The partial wave expansion of the irN scattering amplitude can be worked out straightfor- 
wardly from Eq. ( Appendix A.2 ). By definition, the initial baryon three-momentum p gives 
the positive direction of the z-axis. Inserting the series of Eq. ( Appendix A.2 ) one has for 
the scattering amplitude, 

(7r(-p'; a')N(p', a'; a') \T\n(-p; a)N(p, a; a)} = 
4tt Y e(^)(m(T'(T\£-J)(0aa\£-J)Y £ m (p')T M (s), (Appendix A.4) 

£,m,J 

where T is the T-matrix operator and Tj£ is the partial wave amplitude with total angular 
momentum J and orbital angular momentum £. Notice that in Eq. (Appendix A.4) we 
made use of the fact that V^ m (z) is non-zero only for m = 0. Recall also that because of 
parity conservation partial wave amplitudes with different orbital angular momentum do 
not mix. From Eq. ( Appendix A.4 ) it is straightforward to isolate Tjg with the result, 

1 



Tje(a', a'; a, a 



E / d p' M-p'; a ') Ar (p'' ^; «') l T k(-p; «Mp> «)> (maVi^L)ynp' 



(Appendix A. 5) 
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in the previous expression the resulting Tje is of course independent of choice of a. 
The relation between the Cartesian and charge bases is given by 



|7T + ) = 




+ *k 2 » , 




IO = 




-*k 2 », 




k°> = 


k 3 ) . 




(Appendix A. 6) 



According to the previous definition of states \tt + ) = — 11,+1), |7r~) = |1, — 1) and \tt°) = 
1 7r 3 ) = |1, 0), where the states of the isospin basis are placed to the right of the equal sign. 
Notice the minus sign in the relationship for \ti + }. Then, the amplitudes with well-defined 
isospin, 1 = 3/2 or 1/2, are denoted by Tjji and can be obtained employing the appropriate 



linear combinations of Tjt(a',a';a,a), Eq. ( Appendix A.5 ), in terms of standard Clebsch- 
Gordan coefficients. 

Due to the normalization of the states with well-defined total angular momentum, 
Eq. (Appendix A.3), the partial waves resulting from Eq. (Appendix A.5) with well defined 



isospin satisfy the unitarity relation, 



lmT ue = P \T IM \ 2 . (Appendix A. 7) 

For |p| > and below the inelastic threshold due the one-pion production at |p| ~ 210 MeV. 
Given the previous equation, the S-matrix element with well defined J, J and £, denoted by 
Sjji, corresponds to 

Sue = 1 + i—^-=T Ije , (Appendix A. 8) 

Ally's 

satisfying SueSj je = 1 in the elastic physical region. In the same region we can then write 

S Ije = e 2i5lJl , (Appendix A.9) 

with 5jj£ the corresponding phase shifts. And, form Eqs. ( Appendix A. 8 ) and ( Appendix A.9 ) 
one has 

T IM =^^-smS IM e iSlje . (Appendix A. 10) 

However, if the calculation is perturbative, the S'-matrix does not fulfil unitarity exactly and 



one cannot use Eq. (Appendix A. 10) to calculate the phase shifts. Instead, is necessary to 
perform a perturbative expansion of the previous equation up to the order considered to find 
a relation between the perturbative amplitude to its corresponding phase shift. Following 
this procedure, we find that up to 0(p A ) the different phase shifts can be obtained from the 
perturbative amplitudes by means of the equation: 



UJl 



iPl 



87T^/i 
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ReTjji (Appendix A. 11) 



Appendix B. Tree Level Calculations 

In this Appendix we show the results concerning the tree level calculation of irN scatter- 
ing amplitude. The Born-terms, which are expressed in terms of the Mandelstam variables 
s and u, include also their crossed version. 

Appendix B.l. 0(p) 
• Born-term: 



A + (s,t,u) 
B + (s,t,u) 
A~(s, t, u) 
B~(s, t, u) 



g 2 (m 2 + m N ) 



4/ 2 



s — m 



N 



s — m% 



+ 



u — m 



N 



u — rrin 



4/ 5 



(s + 2m 2 m N + m 2 N ) (u + 2m 2 m N + m 2 N ) 



g 2 (m 2 + m N ) 



[s — rrin 



s — m 



(u - ml) 



4f 



N 



u — m 



N 



s — rrin 



u — m% 



4/ s 



' s + 2m 2 rriN + m 



N) 



(s - ml) 



(u + 2m 2 m N + m 2 N ) 
(u - ml) 



Where m 2 = m — AciM 2 includes the 0(p 2 ) correction to the nucleon mass. 



• Contact term: 



A + (s,t,u) = B + (s,t,u) = A~(s,t,u) = 
B~(s,t,u) = — 

Appendix B.2. 0(p 2 ) 
• Contact term: 



A + (s,t,u) 



B + (s,t,u) 



A (s,t,u) 



B (s,t,u) 



1 

I 2 




-A Cl M 2 + 



im 



+ c 3 (2M 2 -t) 



+ C(p 4 



N 



C4{S — U) 

2c^m^ 

J 77 
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Appendix B.3. 0(p 3 ) 
• Born-term: 



A + (s,t, u) 



B + (s,t,u) 



Ag(2d w - d ls )m N Ml 

J 7T 

Ag(2d w - d 18 )m 2 N Ml(s - u) 



f 2 (s - m%)(u - m 2 N ) 
A~(s,t,u) = 

B _ t = 2g(2d 16 - d 18 )M 2 {3m% - su - m 2 N {s + u)) 
• Contact term: 



(d u - d 15 )(s - u) 
f 2 

A~(s, t, u) = [2(d! + d 2 + 2d 5 )M 2 - {d 1 + d 2 )t + 2d 3 (s - u) 2 ] + 0(g 5 



2m N f 2 
B-(s,t,u) = 



• A Born-term 0(p 3 ^ 2 ): 

An explicit calculation of the Born-term in the s-channel gives: 

A ± = --r^r C t—^r i m N ~ 2 {Ml + 2s) m 3 N - 2m A {Ml + s) 



+ {Ml - AsMl + 3s(s + t)) m N + 2m A {Ml - sf + 3m A st) 
h 2 1 

5 ± = -—^-Cf T (m% - 2 {Ml + 3s) m 2 N - 2m A {m 2 N - Ml + s) m N 

4/~m A s — m A v ' v ' 

+ (M^ - s) 2 + 3sf) (Appendix 
With Cj = 1/9 and Cf = -1/18. 
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A Born-term 0{p 5 ^ 2 ): 



h A ±F A (s,t) 



2/|m A s - m A 
S ± = ^^cf^^ (Appendix B.2) 



2f 2 m\ 



m 



With: 

F A (s, t) = J- (d£ (m% - Ml - s) (m% - 2 (M T 2 + 2s) m% - 2m A (M 2 + s) 

+ (M* - 4sM 2 + 3s(s + t)) m N + 2m A (M 2 - s) 2 + 3m A st) 
+2<#m A + 2 (M 2 + 2s) m% - (M* - 2sM 2 + s(5s + 3t)) 

+m A ((mfy - s) 2 - m^v + s (2 {M 2 - sf + 3st 



m 



N 



F B (s, t) = — (<# (rr& - M 2 - s ) (r< - 2 (M 2 + 3s) m% 
-2m a (m 2 N -Ml + s)m N + {Ml - sf + 3stJ + 2dfm A (m A (3m^ 
-4 (M 2 + s) m 2 N + {Ml - sf + 3si\ - m N (m% - 2 {Ml + 2s) m 2 N 
+M*-4M 2 s + 3s{s + t)))) . 



Appendix C. Loop Level Calculations 

In this section we list the scalar and tensor integrals needed for the one-loop calcula- 
tions performed in this work. These integrals, calculated in dimensional regularization, are 
denoted by T-L mn , where the subscripts m and n correspond to the number of mesonic and 
baryonic propagators, in order, that each integral has. 

We will use the following variables 

w = (p + q y = (p' + q 'Y, 
a^ = { q ' - q y = (p - p'y, 
qh = (p> + py, 

where P (P') corresponds to the incoming (outgoing) nucleon and q (q') to the incoming 
(outgoing) pion, so P 2 = P' 2 = m 2 N and q 2 = q' 2 = M 2 . 
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Appendix C.l. Definitions 
• 1 meson, nucleons: 



i r d d k 



i J (27r) d M 2 — k 2 
M 2 ( 



M 2 ( M 2 

n 10 = 2\M 2 + — -.log — ) (Appendix C.l) 



mesons, 1 nucleon: 



i r d d k 

rtoi 



i 7 (27r) d m 2 -P 

2 / 2 

2 m , / m 



167T 2 \ /I 



Hoi = 2Am i + T7r ^log — ) (Appendix C.2) 



2 mesons, nucleon: 



{"H20) %2o) 



1 /" A 



% J (27r) d (M 2 - k 2 )(M 2 - (k - A) 2 ) 

^20 - -^^2o(t) 

U% = (AW - (?^A 2 )^S(t) + A"A v H^(t) 
1 meson, 1 nucleon: 



Hn(P 2 ) 



-2A 



{1,^} 

(M 2 - k 2 )(m 2 - (P - k) 2 ) 

{ , /u 2 \ P 2 -m 2 + M 2 
1 + log 



16vr 



?7r 



2P 2 



log 



V /4M 2 P 2 - (P 2 - m 2 + M 2 )' 



x 



arctan 



P 2 



m 2 - M 2 - P 2 



arctan 



\ v/4M 2 P 2 - (P 2 - m 2 + M 2 ) 2 

m 2 -M 2 + P 2 
^AM 2 P 2 - (P 2 - m 2 + M 2 ) 2 

H^^P^iP 2 ) 
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AP 



(Appendix C.3) 



• mesons, 2 nucleons: 



i r d d k 

it-02 



% J (2ir) d (m 2 - k 2 ){m 2 — (k — A) 2 ) 
2 mesons, 1 nucleon: 



\n-2i, n-215 ' t 2i J — ~ / 



d d fc {1,^,^} 
(27r) d (M 2 - A; 2 )(M 2 - (k - A) 2 )(m 2 - (P - &) 2 ) 



U% = g^U { £>(t) + Q»Q u U^{t) + A^A^W + ^(AV + Q» A V )H$ (t) 
1 meson, 2 nucleon: 



j-iy -1//' u^X — I 



C^fc {1, fc", Jfe""} 

(2vr) d (M 2 - fc 2 )(m 2 - (P x - A;) 2 )(m 2 - (P 2 - k) 2 ) 



For the topologies displayed in Fig. ® one of the momenta is always on-shell. Choosing 
this momentum to be Pi, we have for the diagram (m): P x = P and P 2 = P'. This case 
defines the integral %A{i) as follows: 

= n 12 (m 2 J t). 

For this case, the tensor decomposition is defined as, 

H%(t) = g^nfit) + Q»Q u nf(t) + AWHffit). 

For the diagrams (c), (d), (g) and (/i) we have instead P\ = P and P 2 = P + q, with the 
integral "Hb(s), 

H b (s) = Hi 2 (s,M 2 ) 
In this case, the tensor decomposition is defined as, 
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• mesons, 3 nucleon: 



H 



03 



d d k 



{27r) d (m 2 - k 2 ){m 2 — (k — Pi) 2 )(m 2 - (k - P 2 ) 2 



This integral can appear in two different configurations. In the first one P± = q and P 2 = q', 
which is labeled in Sees. Appendix C.2 and Appendix C.3 as "Ho3(^ M 2 ). For this case 
its dependence on the t variable comes from the combination (P x — P 2 ) 2 '. For the second 
configuration, however, the t dependence comes from P|, because in this second case P\ = q 
and P2 = q — q' ■ This configuration is labeled in Appendix C.2 and Appendix C.3 as 

n 03 (M 2 ,t). 



1 meson, 3 nucleon: 



{^13,^13} 



t) 



d d k 



i J (27r) d (M 2 - k 2 )(m 2 - (P - A;) 2 )(m 2 - (E - A;) 2 )(m 2 - (P' - fc) 2 ) 



Q»n^(s,t) + (A + 2qyu {2) 



13 
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Appendix C.2. Coefficients of the Passarino-Veltman 
Decomposition 



^ii(s) = f [(a - m 2 + M 2 )^ n (s) + ftio - «oi] 



Uf(t) 
Hf(t) 



2s 

n 11 {m 2 ) + M 2 H A {t)-H 02 {t) 
Am 2 — t 

2M 2 n n (m 2 ) + 2M 2 (M 2 - 4m 2 + t)H A {t) - (2M 2 - 4m 2 + t)H 02 (t) 

2(2 -d) (Am 2 -i) 
M 2 ((l - rf)M 2 + 4m 2 - t)^ A (t) + (1 - d)M 2 U u (m 2 ) 

(2-d)(4m 2 -t) 2 
(2(1 - d)M 2 + (3 - d)(4m 2 - t))U 02 (t) M 2 n u (m 2 ) + ft 10 - U 01 



(Appendix C.4) 



2(2-d)(4m 2 -t) 2 
1 



2(M 4 + (m 2 - s) 2 - 2M 2 (m 2 + s)) 
n n (s)) + (s - m 2 - M 2 )U xl (m 2 ) + 2M 2 U 02 (M 2 ) 



Am 2 (Am 2 - t) 
s-m 2 + M 2 )((s - m 2 - 2M 2 )H B (s) 



+ 



n { 2 \\t) 



2(M 4 + (m 2 - s) 2 - 2M 2 (m 2 + s)) L 

(M 2 - m 2 - 3s)%n(s) + (s + 3m 2 - M 2 )Hu(m 2 

2(s - m 2 )U m (M 2 ) 
2(M 4 + (m 2 - s) 2 - 2M 2 {m 2 + s)) 
(2M 2 - t)U 21 (t) - 2Hn(m 2 ) + 2U 20 (t) 



s - m 2 )(s + 3m 2 - 3M 2 )U B (s) 



2(4m 2 - t) 

-2(M 4 + m 2 (t - AM 2 ))V, 2l (t) + 2(M 2 - 2m 2 )V, ll (m 2 ) + (t - 2M 2 )U 20 (t) 



H%(t) 



2(d - 2) (4m 2 - t) 
(A(d - 1)M 4 - AM 2 ((d - 2)t + 4m 2 ) + t((d - 2)t + 4m 2 ))^ 20 (t) 



+ 



1 



4(c/-2)m 2 (t-4m 2 ) 5 



+ 8m 4 )^ n (m 2 ) + (d- 2) (4m 2 - t)(H i - U 



A(d - 2)(t - Am 2 ) 2 

((d - 2)(M 2 + 2m 2 )t + 4(3 - 2d)M 2 m 2 

(d - l)(t - 2M 2 )V, 2Q (t) 



10 



2(d-2)(t-4m 2 ) 2 
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(AM 2 + 2m 2 



4(M 4 - 2M 2 (s + m 2 ) + m 2 (m 2 - 2s) + s(s + t)) L 
- 2s - t)"Hi 3 (s, t) + 2(m 2 - s - M 2 )n B (s) + (2s + t - 2m 2 - 2M 2 )U A (t) 

+ 



(4M 2 -t)^ 03 (M 2 ,t) 



+ 



4(M 4 - 2M 2 (s + m 2 ) + m 2 (m 2 - 2s) + s(s + t)) 

\l M 2 f2s 4- 1 — 2m 2 ) 

4(M 4 - 2M 2 (m 2 + s) + m 4 - 2m 2 s + s(s + *)) L V 1 j 

(m 2 - s)(4m 2 -t) - 2M 4 )H 13 (s,t) + (2s + t - 2m 2 - 2M 2 )H 03 (M 2 , t) 

2(s + m 2 - M 2 )H B (s) + (t- 4m 2 )U A (t) 
4(M 4 - 2M 2 (m 2 + s) + m 4 - 2m 2 s + s(s + t)) 



Appendix C.3. Results for the loop diagrams 

We list in this section the results concerning to the loop integrals. For the diagrams a + b, 
c + d, e, / , g + h, i, n + o and p + r only the direct version is shown. To construct the full 
contribution of the mentioned diagrams is necessary to add its crossed version according to 
the following rules 34j : 



a $otal( s i t, u) = A ± (s, t) ± A ± (u, t) 
[s,t,u) = B ± (s,t)T B ± (u,t) 



B 



± 

TOTAL 



Loops a+b 



BUS) 



g 2 m 



n i + Hw + (s - m 2 - M 2 )Hn(s) -(s + m 2 )n[\\s 
4(3m 2 + s)(^ i -M 2 Hu{s)) 



8f 4 

B a+ b ( S ) 



rrr 



+ 4(m 2 + s)^ n (s) 
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Loops c+d 



4 

g m 



A t+d( S ) 8/ 4 (s _ m 2; 



[Am 2 n i - 2(m 2 + s)U w 



+ 2(2m 2 (s - m 2 ) + M 2 (m 2 + s))n u (m 2 



- 2(s 5 



m 



M 2 (3m 2 + s))7fn(s) - 4m 2 (s + m 2 )?^™ 



+ (m 4 + 10m 2 s + 5s 2 )ftff (*) + 4m 2 (m 2 (m 2 - 2s) + s 2 )H B (s) 
+8m 2 (m 2 - s)(2(s + m 2 ) - M 2 )] 



,9 J 



[-(s + 7m 2 )H i +4m 2 H 



10 



+ 4m 2 (3m 2 + s)H 02 (M 2 ) + Am 2 M 2 Uu{m 2 ) + (M 2 (s + 3m 2 



+ 4m 2 (s - m 2 ))H n (s) - 8m 4 H { 3(m 2 ) - (s 2 



6sm 2 + ra 4 )7/^(s 



A 



c+c 



- Am 2 M 2 (s + 3m 2 )^ B (s) + 4m 2 (s - m 2 ) 2 wg ) (s) 
+4m 2 (m 2 - s)(3m 2 + s)U ( b\s) 



B- +d {s) 



• Loop e 

A tW = TKfU—2 V [ 2 ( s + 3m2 )^oi + 2( 5 - m 2 )U 10 

lbj^{m z — s) 

+2((s - m 2 ) 2 - M 2 (3m 2 + s))W u (s) - (s - m 2 ) 2 ^!^) 

= 1fi ,4/ 3 ^ 2 \2 [( 9m " + 6m2s + s2 ) n oi + 4m 2 (s - m 2 )Hi 

+(4m 2 (s - m 2 ) 2 - M 2 (3m 2 + s) 2 )7*n(s) - (m 2 - s) 3 ^?^) 

A e -( S ) = A+( S ) 
B-{s) = Bt{s) 
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• Loop f 



m(s — m 2 )'H^(s) 
f W - 2/i 

(m 01 -U 10 -4M 2 U u (s) 



Aj(s) 



A, (a) 



Bj(s) 



2 



8f 4 



Loops g+h 



B g+h (s) 



g 2 m(s — m 2 ) 



2f 



-2H u (s) + Mi 00 + %rn 2 U%\s) 



4f 4 



-2U m + ft 10 + %m 2 H m (M 2 ) + 2M 2 Hn(s) - 2m 2 Mi(m 2 ) 

-2(m 2 + s)Mi } (s) - 8M 2 m 2 H B (s) - 8m 2 (m 2 - s)(7/g } (s) - M^OO) 
= 
= 
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Loop i 

Af(s,t) 



KM 



■^0 [8m 2 (^ 02 (M 2 ) - n 02 (t)) + 2(4m 2 - M 2 )^ n (m 2 ) 

+ (m 2 - s)?4i («) + 2(M 2 - (s + 3m 2 ))H n (s) + 32m 4 (m 2 - s)U 13 (s,t) 
+ 8m 2 M 2 H A (t) - 32m 4 H ( A 1) (t) + 8m 2 (2s + t - 2M 2 - 2m 2 )uf [t) 
- 8m 2 (M 2 + m 2 - s)U B {s) + 8m 2 (M 2 + 3m 2 + s)H%\s) 
+8m 2 (M 2 + m 2 - s)H { b\s) 

1- Hio + 4m 2 (2H 02 (M 2 ) + ft 02 (*)) + 16m 4 ^ 03 (t, M 2 ; 



16f 



+ (M 2 - 4m 2 )-H n (s) - Am 2 wll(m 2 ) + (m 2 - s)^! (*) + 16m 4 M 2 -H 13 (s, t) 
+ 32m 4 (m 2 - s)?^ (M) - 4m 2 M 2 n A (t) + 8m 2 H { J\t) - 8m 2 M 2 H B (s) 



(i) 



Br(M) 



+8m 2 (s + 3m 2 )H i? 

A + (M) 
3 



s) + 8m 2 (m 2 - s)^ 2) (s 



• Loop k 

At(t) = 
Bt(t) = 
^(f) = 

B k W - j 4 



Loop 1 



g 2 m 



V 4 



-2Hoi + 2M 2 Hn(m 2 ) + (Af 2 - 2t) (4m 2 ^ 1 1 ) (t) - ft 20 (t)) 



fl+(t) = 



AT® 



4g 2 m 3 



f 4 
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Loop m 



A+(s,t,u) 
B+(s,t,u) 



A m (s,t,u) 



B m (s,t,u) 








2 ^ 

g m 



s-u)Hf{t) 



9 



f 4 

2 



8f 4 



Hio - Am 2 (H^{m 2 ) - H» 2 {t) + M 2 H A {t) - 2H^{t)) 



(2), 



g 2 m 



f 4 



Loops n+o 

B- +0 (s,t) = B+ +Q (s,t) 



[Hoi - M 2 Hn(m 2 )] 



g 2 (7m 2 + s)(U m - M 2 Hn{m 2 )) 



• Loops p+ 

A p+r (s,t) = 
(M) = 



2^ 

# 2 ( g + 3m 2 )^ 10 
4/ 4 m 2 - s 
Ap +r (s, t) 



• Loops t+u 

At +u {s,t) = 9 ^{H m -M 2 H ll (m 2 )) 
Bt +u {s,t) = ^{U m -M 2 U ll {m 2 )) 

B; +u (s,t) = 
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• Loop v 

A+(s,t) = 
B+(s,t) = 

4T(M) = o 

Appendix D. Identifying the power counting breaking terms 

In this Appendix we explain the method we used to extract analytically the power 
counting breaking terms from the 0(p 3 ) loop amplitude (T) oops ). First, as we did with the 
full amplitude, we decompose Ti oops in terms of its scalar integrals using the Passarino- 
Veltman decomposition. 

T[oops ^ ^ ^mnH-mni 



where the scalar integrals T-L mn are defined in |Appendix C| and C mn refers to its coefficients 



that result in the Passarino-Veltman decomposition. Second, we calculate the infrared regu- 
lar part (26[ of these scalar integrals (TZ mn ), because it contains all the PCBT. Its calculation, 
for each of the scalar integrals, is straightforward because the chiral expansion of the regular 



part commutes with the integration in the Feynman parameters [85j. The chiral order of 
each (suin tell us up to which order in the chiral expansion we need to obtain the regular 
part of H m n- So, finally, we obtain, 

Tioops = y^^mnUmn- (Appendix D.l) 



We expand Eq. ( Appendix D.l ) in a chiral series to end with a string of terms that can 
be splitted into a part that has chiral order lower than three (these are the PCBT in our 
case) and an infinite series that respect the power counting. The (finite) terms that break 
the power counting have the same analytical structure than the monomials in the original 
Lagrangian and can be cancelled via a LECs redefinition (see the next Appendix). 



Appendix E. Low-Energy Constants Renormalization 

In this appendix we show how to redefine the 0(p 2 ) and 0(p 3 ) LECs in order to cancel 
the divergences and the PCBT. In this way we have full relativistic scale-independent chiral 
amplitudes free from divergences that respect the chiral power counting. 
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Appendix E.l. 0(p 2 ) LECs 

The 0(p 2 ) LECs are redefined in order to cancel both divergent parts, as well as PCBT. 



m 
~7 



Ci —> c 

C2 _ 40MS + 2 xil^pl + _5L_ [(9 » - i,» log g ) - (2 + 9 4 )] 



(# 4 -6# 2 + l)m 



c 3 + 2A - - 2 ' ^ + ^[(<7 4 - 6, 2 + 1) log - ) + V 



c 4 - cf — + 2A (3 ^ y 1)m + ^[(3, 4 - 2/ - 1) log g'] - , 2 (5 + , 2 )] 



Appendix E.2. 0(p 3 ) LECs 

In contrast to the C(p 2 ) LECs, the C(p 3 ) ones only cancel divergent parts (along with 
their scale-dependent logarithms) because the 0(p 3 ) analytical terms do not break the power 
counting in our 0(p 3 ) calculation. They are renormalized within the MS scheme (also known 
as WS — 1). 

, xm^ T 3# 4 -V + l 3£ 4 -V + 1 (m 2 
d l + d 2 --r (dt + d 2 ) MS + 2A y + \ ao y — log ' 



48/ 2 768vr 2 / 2 °\fi 2 

1 + 8 g 2 + 8 < 
d 5 df b - 2A— — - — 9J9 log 



ims n \9 +« 5+8 /m' 
/s " 2A l87^~76W 

di4 - d 1B -> (du - rfisj + — + 647r 2 y 2 io S I ^5" 
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Appendix F. Summary of results for A-ChPT in the IR scheme 



In Table IFTT71 we show the results at \fs max = 1-20 GeV of the LECs and the Xd.o.f. f° r 
A-ChPT within the IR scheme, which is to be compared with the EOMS results in Table [2J 
In Table IF. 181 we give the results obtained at this energy for a selection of observables. For 
the calculation of Act we have used the methods shown in Sec. 15. 2\ whereas for the rest 
of the observables we have made use of the formulas for d~Q , and ovjv given for IR in 



Refs. 26, 34 



LEC 


KA85 A-IR 


WI08 A-IR 


ci 


-0.196(31) 


-0.371(30) 


C2 


1.88(10) 


1.97(9) 


C3 


-2.90(13) 


-3.16(12) 


C 4 


1.81(6) 


1.96(6) 


di + d 2 


0.91(9) 


1.23(9) 


ds 


-1.26(6) 


-1.36(6) 


d 5 


0.168(42) 


-0.124(40) 


du — d±5 


1.33(16) 


-0.97(15) 


dm 


-2.66(26) 


-1.79(25) 


h A 


3.096(35) 


2.956(34) 


Xd.o.f. 


4.15 


1.67 



Table F.17: Result for the LECs obtained in a fit up to energies of \fs n 
A-ChPT. 



1.20 GeV in IR 





KA85 A-IR 


WI08 A-IR 




8.1(8)% 


5.5(8)% 




-0.9(12) 


-0.81(11) 




0.338(41) 


0.420(38) 


a wN (MeV) 


-5.8(2.4) 


7.7(2.3) 



Table F.18: Summary of the results obtained in A-IR for some of the observables studied in this paper. 



61 



References 



[1] G. Hoehler, in Landolt-Bornstcin, 9b2, cd. H. Schoppcr (Springer, Berlin, 1983) 

[2] J. Gasser, M. E. Sainio and A. Svarc, Nucl. Phys. B 307 779 (1988). 

[3] E. Epelbaum and U. -G. Meissner, larXiv:1201."2l36l [nucl-th]. 

[4] Computer code SAID, online program at |http://gwdac.phys.gwu.edu/~| solution WI08. R. A. Arndt 

et at, Phys. Rev. C 74 045205 (2006). solution SM01. 
[5] E. Matsinos, W. S. Woolcock, G. C. Oades, G. Rasche, A. Gashi, Nucl. Phys. A 778 95 (2006). 
[6] V. Baru, C. Hanhart, M. Hoferichter, B. Kubis, A. Nogga, D. R. Phillips, Phys. Lett. B 694 473 

(2011). 

[7] M. C. M. Rentmeestcr, R. G. E. Timmermans, J. L. Friar and J. J. de Swart, Phys. Rev. Lett. 82 
4992 (1999);M. C. M. Rentmeestcr, R. G. E. Timmermans and J. J. dc Swart, Phys. Rev. C 67 044001 
(2003); D. R. Entcm and R. Machlcidt, Phys. Rev. C 68 041001 (R) (2003); E. Epelbaum, W. Glockle 
and U.-G. Mcifincr, Nucl. Phys. A 747 362 (2005); Eur. Phys. J. A 19 125 (2004). 
[8] S. R. Beane, W. Detmold, K. Orginos and M. J. Savage, Prog. Part. Nucl. Phys. 66, 1 (2011). 
[9] M. Fukugita, Y. Kuramashi, M. Okawa, H. Mino and A. Ukawa, Phys. Rev. D 52, 3003 (1995). 
[10] S. R. Beane, W. Detmold, T. C. Luu, K. Orginos, A. Parreno, M. J. Savage, A. Torok and A. Walker- 
Loud, Phys. Rev. D 79, 114502 (2009). 
[11] S. Weinberg, Physica A 96 327 (1979). 
[12] J. Gasser and H. Leutwyler, Annals Phys. 158 142 (1984). 

[13] J. Gasser and H. Leutwyler, Annals Phys. 158, 142 (1984); Nucl. Phys. B 250, 465 (1985). 
[14] H. Leutwyler, Annals Phys. 235, 165 (1994). 

[15] G. Eckcr, Prog. Part. Nucl. Phys. 35, 1 (1995); A. Pich, Rept. Prog. Phys. 58, 563 (1995); S. Scherer, 
Adv. Nucl. Phys. 27, 277 (2003); U.-G. Mcifincr, Rep. Prog. Phys. 56 903 (1993). 

[16] Dashen R.F. Adlcr, S.L. Current Algebras and Applications to Particle Physics. New York, USA: 
Benjamin 394 p (1968). 

[17] M. L. Goldberger and S. B. Trciman, Phys. Rev. 110, 1178 (1958); M. L. Goldberger and S. B. Trciman, 

Phys. Rev. Ill, 354 (1958). 
[18] T. P. Cheng, R. F. Dashen, Phys. Rev. Lett. 26 594 (1971) . 
[19] S. Weinberg, Phys. Rev. Lett. 18, 188 (1967). 

[20] S. L. Adlcr, Phys. Rev. 137 (1965) B1022. ibid 139 B1638 (1965). 
[21] V. Bernard, N. Kaiser and U.-G. Mcifincr, Int. J. Mod. Phys. E 4 193 (1995). 
[22] V. Bernard, Prog. Part. Nucl. Phys. 60, 82-160 (2008). 
[23] E. E. Jenkins and A. V. Manohar, Phys. Lett. B 255 558 (1991). 
[24] N. Fettes, U.-G. Meifiner and S. Steininger, Nucl. Phys. A 640 199 (1998). 
[25] N. Fettes and U.-G. Meifiner, Nucl. Phys. A 676 311 (2000). 
[26] T. Becher and H. Leutwyler, Eur. Phys. J. C 9 643 (1999). 
[27] T. Fuchs, J. Gegelia and S. Scherer, Eur. Phys. J. A 19, 35 (2004). 
[28] B. R. Holstein, V. Pascalutsa and M. Vanderhaeghen, Phys. Rev. D 72, 094014 (2005). 
[29] L. S. Ceng, J. Martin Camalich, L. Alvarez-Ruso and M. J. Vicente- Vacas, Phys. Rev. Lett. 101 222002 
(2008). 

[30] J. Martin Camalich, L. S. Gcng, M. J. Vicente Vacas, Phys. Rev. D82, 074504 (2010). 
[31] P. J. Ellis and H. B. Tang, Phys. Rev. C 57 3356 (1998). 
[32] K. Torikoshi and P. J. Ellis, Phys. Rev. C 67 015208 (2003). 

[33] J. M. Alarcon, J. M. Camalich, J. A. Oiler and L. Alvarez-Ruso, Phys. Rev. C 83 055205 (2011); we 
have found numerical errors in the results published in this paper. An erratum and a revised version 
of the pre-print are in due course. 

[34] T. Becher and H. Leutwyler, JHEP 0106 (2001) 017. 

[35] T. Ledwig, J. Martin-Camalich, V. Pascalutsa and M. Vanderhaeghen, Phys. Rev. D 85, 034013 (2012). 
[36] V. Pascalutsa. liarXiv: 11 10.57921 [nucl-th], 

[37] E. E. Jenkins and A. V. Manohar, Phys. Lett. B 259, 353 (1991). 

[38] T. R. Hemmert, B. R. Holstein, J. Kambor, J. Phys. G G24, 1831-1859 (1998). 



62 



[39] V. Pascalutsa, M. Vanderhaeghen and S. N. Yang, Phys. Rept. 437 125 (2007). 

[40] N. Fcttes and U.-G. Meifiner, Nucl. Phys. A 679 629 (2001). 

[41] J. Gcgelia and G. Japaridzc, Phys. Rev. D 60, 114038 (1999) . 

[42] T. Fuchs, J. Gegelia, G. Japaridzc and S. Schcrcr, Phys. Rev. D 68, 056005 (2003). 

[43] V. Pascalutsa and D. R. Phillips, Phys. Rev. C 67, 055202 (2003). 

[44] R. Koch, Nucl. Phys. A 448 (1986) 707; R. Koch and E. Pictarincn, Nucl. Phys. A 336 331 (1980). 

[45] S. Weinberg, Phys. Lett. B 251 288 (1990); Nucl. Phys. B 363 3 (1991). 

[46] J. A. Oiler, M. Verbeni and J. Prades, JHEP 0609 079 (2006). 

[47] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321, 311 (1989). 

[48] V. Bernard, N. Kaiser and U.-G. Meifiner, Nucl. Phys. A 615 483 (1997). 

[49] V. Pascalutsa, Phys. Rev. D 58, 096002 (1998); V. Pascalutsa and R. Timmermans, Phys. Rev. C 60, 

042201 (1999); V. Pascalutsa, Phys. Lett. B 503, 85 (2001). 

[50] H. Krcbs, E. Epelbaum and U. -G. Meissner, Phys. Rev. C 80, 028201 (2009). 

[51] L. S. Ceng, J. Martin Camalich, L. Alvarez-Ruso, M. J. Vicente Vacas, Phys. Rev. D78, 014011 (2008). 

[52] K. Nakamura et al. [Particle Data Group Collaboration], J. Phys. G G 37 075021 (2010). 

[53] S. Weinberg, Cambridge, UK: Univ. Pr. (1996) 489 p. 

[54] N. Fcttes and U. -G. Meissner, Nucl. Phys. A 693, 693 (2001). 

[55] V. Pascalutsa, M. Vanderhaeghen, Phys. Lett. B 636 31 (2006). 

[56] S. -i. Ando and H. W. Fearing, Phys. Rev. D 75, 014025 (2007). 

[57] M. R. Schindler, T. Fuchs, J. Gegelia and S. Scherer, Phys. Rev. C 75, 025202 (2007). 

[58] T. Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118, 153 (1999). 

[59] B. Long and V. Lensky, Phys. Rev. C 83 045206 (2011). 

[60] B. Tromborg, S. Waldenstrom and I. Ovcrbo, Phys. Rev. D 15, 725 (1977). 

[61] A. Gashi, E. Matsinos, G. C. Oades, G. Rasche and W. S. Woolcock, Nucl. Phys. A 686, 447 (2001). 

Nucl. Phys. A 686, 463 (2001). 

[62] K. H. Augcnstcin, G. Hohler, E. Pietarinen and H. M. Staudcnmaier, Karlsruhe Data Tape, ZAED 

Physics Data (1977). 

[63] The full data base can be accessed via the website |http://gwdac.phys.gwu.edu 

[64] E. Matsinos and G. Rasche. larXiv:1203.3625l 

[65] P. Buettiker and U.-G. Meifiner, Nucl. Phys. A 668 97 (2000). 

[66] J. A. Oiler, E. Oset, J. R. Pelaez, Phys. Rev. D 59 074001 (1999). 

[67] J. A. Oiler and E. Oset, Phys. Rev. D 60, 074023 (1999). 

[68] U.-G. Meifiner and J. A. Oiler, Nucl. Phys. A 673 311 (2000). 

[69] J. A. Oiler and U.-G. Meifiner, Phys. Lett. B 500 263 (2001). 

[70] J. M. Alarcon, J. Martin Camalich, J. A. Oiler, Phys. Rev. D 85, 051503(R) (2012). 

[71] J. M. Alarcon, PhD Thesis (2012). University of Murcia. 

[72] G. C. Oades, G. Rasche, W. S. Woolcock, E. Matsinos and A. Gashi, Nucl. Phys. A 794, 73 (2007) . 

[73] M. Hoferichter, B. Kubis and U. -G. Meissner, Phys. Lett. B 678, 65 (2009); M. Hofcrichter, B. Kubis 

and U. -G. Meissner, Nucl. Phys. A 833, 18 (2010). 

[74] S. Weinberg, Phys. Rev. Lett. 17, 616 (1966). 

[75] V. Bernard, N. Kaiser and U.-G. Meifiner, Phys. Lett. B 309, 421(1993). 

[76] J. J. de Swart, M. C. M. Rentmeester and R. G. E. Timmermans, nN Newsletter 13 96 (1997). 

[77] I. Jameson, A. W. Thomas and G. Chanfray, J. Phys. G G 18, L159 (1992). 

[78] J. Gasser, H. Lcutwylcr, M. E. Sainio, Phys. Lett. B 253, 252 (1991). Phys. Lett. B 253, 260 (1991). 

[79] M. M. Pavan, I. I. Strakovsky, R. L. Workman, R. A. Arndt, PiN Newslett. 16 110 (2002). 

[80] M. G. Olsson, Phys. Lett. B 482, 50 (2000). 

[81] J. M. Alarcon, L. S. Geng, J. Martin Camalich, J. A. Oiler. larXiv:1209.2870l [hep-ph]. 

[82] C. Ditsche, M. Hofcrichter, B. Kubis and U. -G. Meissner, JHEP 1206, 043 (2012) 

[83] V. Bernard, N. Kaiser, U.-G. Meifiner, Phys. Lett. B 389 144 (1996). 

[84] J. A. Oiler and L. Roca, Phys. Lett. B 651, 139 (2007); J. F. Donoghue, J. Gasser and H. Leutwyler, 
Nucl. Phys. B 343, 341 (1990). 



63 



[85] J. Gegelia, G. S. Japaridze and K. S. Turashvili, Theor. Math. Phys. 101, 1313 (1994). 



64 



